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Within the context of supersymmetric space-time (D-particle) foam in string/brane-theory, we 
discuss a Finsler-induced Cosmology and its implications for (thermal) Dark Matter abundances. 
This constitutes a truly microscopic model of dynamical space-time, where Finsler geometries arise 
naturally. The D-particle foam model involves point-like brane defects (D-particles), which provide 
the topologically non-trivial foamy structures of space-time. The D-particles can capture and emit 
stringy matter and this leads to a recoil of D-particles. It is indicated how one effect of such a 
recoil of D-particles is a back-reaction on the space-time metric of Finsler type which is stochastic. 
We show that such a type of stochastic space-time foam can lead to acceptable cosmologies at late 
epochs of the Universe, due to the non-trivial properties of the supersymmetric (BPS like) D-particle 
^ <— j defects, which are such so as not to affect significantly the Hubble expansion. The restrictions placed 

on the free parameters of the Finsler type metric are obtained from solving the Boltzmann equation 
in this background for relic abundances of a Lightest Supersymmetric Particle (LSP) dark matter 
candidate. It is demonstrated that the D-foam acts as a source for particle production in the 
Boltzmann equation, thereby leading to enhanced thermal LSP relic abundances relative to those in 
the Standard ACDM Cosmology. For D-particle masses of order TeV, such effects may be relevant 
for dark matter searches at colliders. The latter constraints complement those coming from high 
■ energy gamma-ray astronomy on the induced vacuum refractive index that D-foam models entail. 

^ 1 We also comment briefly on the production mechanisms of such TeV-mass stringy defects at colliders, 

which, in view of the current LHC experimental searches, will impose further constraints on their 
Q\ ■ couplings. 

o 

(Sj ■ I. INTRODUCTION 

There has been recent interest in physics beyond the Standard Model which can incorporate Lorentz symmetry 
violations |l| in a Finsler geometry setting 0tJ|- The Finsler metrics are functions not only of the space-time 
coordinates but also of the tangent vectors (momenta) at points of the curved manifolds. Physical examples of Finsler 
metrics have readily arisen in analogue condensed matter systems [H and so recent attempts to base Cosmology on 
Finsler have been necessarily purely phenomenological Q. Hence the fundamental viability, from a physical point of 
' view, cannot be tested in such studies. 

In a recent Letter Q we have proposed a scheme, which involved back reaction of space-time defects within a 
supersymmetric string/brane framework. The back reaction has the form of an induced Finsler metric Q, but of a 
stochastic kind. In this longer article, we shall give details and address the implications of the induced stochastic 
Finsler geometry on the Hubble expansion of the Universe and the Dark Matter (thermal) relic abundance. Just as 
candidates for Dark Matter have to be sought from Physics beyond the Standard Model, we need to have a framework 
which goes beyond local effective Lagrangians, to incorporate the Finsler dynamical background. Furthermore, the 
localised space-time defects (D-particles [8-11]) that we consider, are generic in string/brane theory. These localised 
defects, although point-like, nevertheless cannot be just considered as another form of particle excitation, as we shall 
discuss below. We give arguments which show that, owing to the presence of D-particles and D-particle/string bound 
states in the vicinity of the brane Universe, a viable and plausible Cosmology, with small but measurable in principle 
implications for the Dark matter relic abundances, exists. Together with the results of our earlier Letter [7], this 
represents the first attempt to construct a microscopic Finsler Cosmology, based on our current understanding of 
Gravitational Physics beyond the Standard Model, and in particular within a string/brane framework. 

The nature of the Dark sector of our Universe constitutes one of the major unresolved puzzles of modern physics. 
Indeed, according to observations over the past twelve years, 96% of our Universe energy budget consists of unknown 
entities: 23% is Dark Matter (DM) and 73% is Dark Energy (DE), a mysterious form of ground state energy. DE 
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is believed to be responsible for the current-era acceleration of the Universe. These numbers have been obtained by 
best- fit analyses of a plethora of astrophysical data to the so-called Standard Cosmological Model (ACDM), which 
is a Friedmann-Robertson- Walker (FEW) cosmology, involving cold DM as the dominant DM species, and a positive 
cosmological constant A > 0; the data range from direct observations of the Universe acceleration, using type- 
Ia supernovae |12], to cosmic microwave background (l3l. Il4j|. baryon oscillation [l5| and weak lensing data [l6j . It 
should be stressed that the afore-mentioned energy budget depends crucially on the theoretical model for the Universe 
considered. 

An interesting, and not commonly discussed, class of cosmological models that may lead to modifications of the Dark 
sector, involves space-time with a "foamy" structure at microscopic (Planckian or string) scales pjj, due to quantum 
gravitational interactions. In the past, for a variety of reasons, such models (differing in the details of the constructs 
of space-time foam) have been considered by many authors. They exhibit a profusion of features that can be falsified 
experimentally and their predictions range from light-cone fluctuations caused by stochastic metric fluctuations [l8j . 
to macroscopic Lorentz symmetry violations As stated earlier, in this note we would like to present a first study 

examining the contribution of such space-time metric stochastic fluctuations on the Dark sector and in particular on 
the DM sector of the Universe. We shall focus on a particular class of stochastic space-time foam models, inspired 
by certain types of string theory. These involve localized space-time defects (DO-branes or D-particles) [8] which are 
either allowed background configurations in (supersymmetric) type IA string models (8l-[l0| or arise effectively from 
suitable compactifications of higher dimensional branes (e.g. D3-branes wrapped up in appropriate three cycles in the 
context of type IIB strings [llj). Observers on the brane detect a foam-like structure due to the crossing of the brane 
by D-particles. In this higher-dimensional geometry, only gravitational fields are allowed to propagate in the bulk; 
all other particle excitations, including DM candidates, are assumed to be described by open strings with their ends 
attached to the brane world. The brane is assumed to have three large spatial dimensions, and, depending on the 
model of string theory considered, it may have a number of compactified extra dimensions. 

Dynamical D-particles should not be viewed as material excitations of the vacuum but rather as vacuum struc- 
tures. This contrasts with attempts to represent such D-particles as ordinary superheavy DM excitations from the 
vacuum poj . owing to the completely localized nature of the D-particles in the extra dimensions. In our construction, 
closest in spirit to weak coupling string theory, they are just vacuum defects. In fact it can be shown that the gravi- 
tational interactions among such (BPS) D-particles are cancelled by appropriate gauged repulsive forces induced on 
them from other branes in our supersymmetric models of D-foam [8, 9]. Hence such a collection of D-particles in the 
bulk does not affect the Hubble expansion on the D3-brane worlds, and so their concentration cannot be restricted 
by considerations on overclosure of the Universe within our models. However important constraints on the density 
of defects in D-foam models can still be imposed by astrophysical experiments on the arrival times of high energy 
cosmic photons (2l| . According to these string- foam models, there should be an effective refractive index in vacuo, 
such that higher ener gy p hotons would be delayed more, since they would cause stronger disturbance on the back- 
ground space-time pi. I ldL [Tl|. The interaction of material open strings (such as photons) with the D-particles leads 
to an induced distortion of space-time described by a metric, which depends on both the coordinate and momentum 
transfer of the photon during its scattering with the defect and so has similarities to a Finsler metric 0|. This is a 
topologically non-trivial process, involving the creation of a non-local intermediate string state, oscillating from zero 
length to a maximum one, according to a time-space stringy uncertainty [ToL H3|. This causes a time delay for the 
photon emerging after capture by the defect, proportional to the incident energy of the photon. 

The purpose of this article is to analyze, in the same spirit, the modification of the estimate of the DM budget of 
the Universe as compared to the ACDM model due to the quantum fluctuations of the D-particles. 1 In fact, as we 
shall argue, the propagation of massive DM particles on such space-times, induces a back-reaction, which in turn has 
consequences for the amount of thermal DM relics of these particles; this, in turn, impacts on astroparticle tests of 
particle physics models incorporating supersymmetry (SUSY), that provide currently one of the leading DM candidate 
species, the neutralino. Its thermal abundance, calculated within the simplest SUSY models (minimal supersymmetric 
model embedded in minimal supergravity [23]), is heavily restricted by cosmic microwave background data; hence the 
available parameter space for these simplest supersymmetric models may vanish in the near future on incorporating 
also data from collider experiments such as the LHC at CERN [2J|. These constraints depend strongly on details 



The D-foam contributions to the Dark Energy budget have been discussed already in J^, I , and found to be phenomenologically 
acceptable in magnitude, and almost constant for the late eras of the Universe; the contribution may be guaranteed by tuning the 
adiabatic velocity of the slowly moving D-brane world in the bulk. The existing supersymmetries in the construction of [£J] guarantee 
a vanishing vacuum energy when no relative motion between D-branes and D-particles is present. This non-trivial input to DE comes 
from stretched strings between the bulk D-particles and the brane world, in directions perpendicular to the (uncompactified) brane 
dimensions. Any motion of D-particles parallel to the brane world, including their recoil during interaction with matter/radiation on 
the brane, does not lead to vacuum energy contributions. 
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of theoretical models. In the presence, for instance, of extended supergravity models with time-dependent dilaton-</> 
sources [25j . which characterize certain (non-equilibrium) string cosmology models [26], the calculated amount of 
thermal neutralino relic abundance can be smaller than the one calculated within the ACDM-minimal supergravity 
cosmology. Such dilaton models allow more scope for SUSY, which can thus survive the otherwise stringent tests at 
the LHC [27(. In our work we will find that the effects of the D-foam on the thermal relic abundances oppose those 
from the dilaton models. These effects become bigger for models allowing low string mass scales. 

For clarity we will first outline the basic reason behind such modifications in the DM thermal abundances [7]. Non- 
equilibrium cosmology models are associated with space-time distortions, due to either the pre sence of time-dependent 
dilaton sources (cf. string cosmologies involving supercritical (SSC) dilaton quintessence |25|), or the induced back- 
reaction of the DM particles onto the space-time itself. Boltzmann equations are used to determine the thermal 
DM species cosmic abundances. The effect of space-time distortions appears as extra contributions to a source T(t) 
for particle production (at cosmic Robertson- Walker time t) on the right-hand side of the appropriate Boltzmann 
equation. In general (for a Universe with three spatial dimensions (i.e. a 3-brane), with no disappearance of particles 
into the bulk), the apposite Boltzmann equation for n(t), the density of the DM species, reads: 

^ +3Hn = T(t)n + C[n], (1.1) 

da 

where H = — = - is the Hubble ratio, a being the scale factor in the Robertson- Walker metric, and C\n] is the 

a a ' ° i L J 

standard collision term describing the deviations from thermal equilibrium of the species population. In keeping with 
standard cosmologies, this has the form [28] 

C[n] = -(av){n 2 -n 2 cq ), (1.2) 

with n cq a thermal equilibrium density of a heavy dark matter species, a the total cross-section evaluated in the 
background metric, and v the Mueller velocity. 

The cosmological implications of the D-particle foam model are not restricted to Dark Matter abundances. It is 
also natural to explore the possible role of D-particle foam [H, 0] on particle statistics, extending and completing 
the ideas outlined in [29j. The implication in the early universe of the statistical description of particles has been 
investigated [13] within the context of non-extensive statistics pioneered by Tsallis [3l| . It is interesting, owing to the 
pervasive nature of space-time foam within our current physical frameworks, to compare its implications for statistics 
with the paradigm of non-extensive statistics. 

The structure of the article is as follows: in section [TT] we review the model of stringy space-time foam (D-particle) 
on which we base our analysis. We pay particular attention in discussing the induced modifications to the space-time 
metric as a result of the interactions of matter (and radiation) particle excitations with the D-particle defects of 
the foam. As already mentioned, such metrics depend on the momentum variables of the particle, in addition to 
any coordinate dependence, and hence the induced metric is of Finsler type [2]. We also discuss briefly non-trivial 
issues associated with the breaking or better obstruction of target-space supersymmetry that should characterize 
any phenomenologically realistic model. In section ITTT1 we analyze the implications of D-foam on particle statistics. 
The reasons why the statistical distributions of particles change in the presence of the foam is the modification of 
the energy dispersion relations of particle probes, as a result of the stochastic metric background. This affects the 
relevant distributions. An interesting feature of our model is that the foam effects are different between fermion 
and boson excitations. In section ITVl we discuss cosmological aspects of the D-foam, and in particular its effects on 
the Hubble expansion of the Universe and the induced modifications to the Boltzmann equation that determines the 
(thermal) Dark Matter relic abundances. The modifications to the Boltzmann equation are due to the Finsler 0] 
nature of the foam-induced space-time metric. We find that the presence of the foam acts as a source for particle 
production and affects the form of the Boltzmann equation in such a way as to lead to an increase of the DM thermal 
relic abundance in comparison to the foam-free Standard Cosmology model. Moreover, contrary to conventional 
string/brane approaches, we argue that the D-particles in the models of D-foam discussed here behave as a dark 
energy fluid rather than dark matter. Some elementary phenomenological considerations, concerning falsification 
of foam models involving low-string scales via WMAP constraints, are presented in section [V] We also comment 
briefly on the production mechanisms of TeV- mass-scale D-particles at colliders, and argue that current experimental 
searches at the LHC can impose further constraints on the effective couplings of these defects to Standard Model 
particles. The strength of such couplings, however, as well as the production mechanism itself, are highly string- 
model-dependent, and thus we reserve a detailed analysis along these lines for a future work. Finally, in section IVll 
we give our conclusions and outlook. Technical aspects of the work are presented in two appendices. 
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II. STOCHASTIC D-FOAM BASICS 



A class of model stochastic space-time foam has been suggested which is based on a gravitational foam consisting of 
real (rather than virtual) space-time defects in higher dimensional space-time, consistent with the viewpoint that our 
world is a brane hyperspace embedded in a bulk space-time [H, In general the construction of a model involves a 
number of parallel brane worlds with three large spatial dimensions, the required number being determined by target 
space SUSY. These brane worlds move in a bulk space-time containing a gas of point-like bulk branes, called D- 
particles, which are stringy space-time solitonic defects (cf. fig. [TJ . One of these branes is the observable universe. On 
this brane the D-particles will appear as space-time defects. Typically open strings interact with D-particles and satisfy 
Dirichlet boundary conditions when attached to them. Closed and open strings may be cut by D-particles. D-particles 
are allowed in certain string theories such as bosonic and type IIA. Here we will consider them to be present in string 
theories of phenomenological interest [l(| since, even when elementary D-particles cannot exist consistently, as is the 
case of type IIB string models, there can be effective D-particles formed by the compactification of higher dimensional 
D-branes [11] . Moreover D-particles are non-perturbative constructions since their masses are inversely proportional 
to the the string coupling g s . The study of D-brane dynamics has been made possible by Polchinski's realization 
that such solitonic string backgrounds can be described in a conformally invariant way in terms of world sheets with 
boundaries (3^. On these boundaries Dirichlet boundary conditions for the collective target-space coordinates of 
the soliton are imposed [33]. When low energy matter given by a closed string propagating in a (d+ l)-dimensional 
space-time collides with a very massive D-particle (DO-brane) embedded in this space-time, the D-particle recoils as 
a result [3^ . 




Figure 1: Schematic representation of a generic D-particle space-time foam model. The model of ref. [9(, which acts 

as a prototype of a D-foam, involves two stacks of D8-branes, each stack being attached to an orientifold plane. 
Owing to their special reflective properties, the latter provide a natural compactification of the bulk dimension. The 
bulk is punctured by DO-branes (D-particles), which are allowed in the type IA string theory of The presence of 
a D-brane is essential due to gauge flux conservation, since an isolated D-particle cannot exist. Open strings live on 
the brane world, representing Standard Model Matter and they can interact in a topologically non-trivial way with 
the D-particle defects in the foam. Recoil of the D-particle during such interactions creates appropriate distortion in 
the space-time geometry, which depend on the momenta of the incident string states, and thus are of a generalized 

Finsler type. 

In this article, we shall consider the simple case of bosonic stringy matter coupling to D-particles and so we can only 
discuss matters of principle and ignore issues of stability. However we should note that an open string model needs to 
incorporate for completeness, higher dimensional D-branes such as the D3-brane. This is due to the vectorial charge 
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carried by the string owing to the Kalb-Ramond field. Higher dimensional D-branes (unlike D-particles) can carry the 
charge from the endpoints of open strings that are attached to them. For a closed bosonic string model the inclusion 
of such D-branes is not imperative. The requirement of higher dimensional branes is for consistency, but is, otherwise, 
not pertinent to our analysis. The current state of phenomenological modelling of the interactions of D-particle foam 
with stringy matter will be briefly summarized now. Since there are no rigid bodies in general relativity the recoil 
fluctuations of the brane and their effective stochastic back-reaction on space-time cannot be neglected. D-particle 
recoil in the "tree approximation", i.e. in lowest order in the string coupling g s , corresponds to the punctured disc or 
Riemann sphere approximation in open or closed string theory and induces a non-trivial space-time metric. 



A. Elements of Logarithmic Conformal Field Theory of D-particle Recoil and Space-Time Implications 

In a flat target-space time, the recoil conformal field theory has been developed in [34],[35||. The pertinent vertex 
operator for our purposes reads: 

v= f u i x°e e (x°)d n x i (2.1) 

with 

e*(*) = r -rt^. (2.2) 

Above, Ui is the recoil velocity of the D-particle, related to the momenta (fc( 2 )) of the propagating string state 
before (after) the recoil via: 

Ui=S< ^ = ^ (2 - 3) 

as a result of momentum conservation during the scattering event. In the above relation Afej is the momentum transfer 
during a collision, g s < 1 is the (weak) string coupling, and M s is a string mass scale. M& = M s /g s is the D-particle 
mass, which may be identified with the Planck mass Mp. These relations have been calculated for non-relativistic 
DO-branes where Uj is small. 

The operator (|2.ip has non-trivial anomalous dimension — e 2 /2 < 0, and as such the theory requires Liouville 
dressing to restore the conformal symmetry It can be shown that the theory is supercritical (once the underlying 
string theory lives in its critical target-space dimension) , and the Liouville dressing results in the operator (|2.1|) being 
replaced by the world-sheet bulk operator 
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a(p Ul (e e (t) d a td a X i + t5 E (t)d a td a t + td 2 X i ) , (2.4) 



where a = is the Liouville anomalous dimension. Upon considering times t relatively long after the scattering, 
and using the world-sheet equations of motion in flat target space time, we observe that only the first term remains. 
Moreover, upon identifying the target-time t with the Liouville zero mode ipo/y/2, we obtain for the target-space time 
metric pi l3rj| 

9ij = <%> 3oo = -1, goi = ^Uie^~ et = ~Ui , t > 0, (2.5) 

where the suffix denotes temporal (Liouville) components. Notice that the world-sheet cr-model time t is assumed 
of Euclidean signature, for reasons of convergence of the world-sheet path integral. It is in this signature that 
the logarithmic conformal algebra of recoil is rigorously valid [34], However, the target-space time signature is 
Minkowskian, due to the supercriticality of the Liouville string, after the identification of the Liouville mode with the 
target time (ip — s/2t): 

ds 2 = -dip 2 + dt 2 + u t dtdX l + (dX J ) 2 = -dt 2 + udtdX 1 + {dX J ) 2 . (2.6) 



The fact that the induced deformation goi depends on the energy content (ui, cf. (|2.3|l ) of the low energy particle, 
implies [8|] that the space-time metric in this case is akin to a Finsler metric [2| . Such a feature does not arise in other 
approaches to space-time foam and will be important in the formulation of our microscopic model. In a stochastic 
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D-particle recoil model the velocity m is taken to be random (owing to a lack of any particular prior knowledge). 
Specifically it is taken (38| to have a distribution which is Gaussian with zero mean and variance a. The zero mean 
of the distribution preserves Lorentz symmetry on average, which is a desirable feature for a (string theory) model, 
both conceptually and phenomenologically. However, the stochastic recoil velocity fluctuations break the symmetry 
(but the effects are suppressed since the respective variances are assumed to be inversely proportional to the effective 
Planck mass) and, on cosmological scales, macroscopic effects of such a breaking may arise. The form of this stochastic 
induced metric will serve as the starting point of our approach to non-conventional statistics [29]. 

Embedding the space-time D-foam in a Friedmann- Robertson- Walker expanding Universe with scale factor a(t), 
requires some discussion. First of all we need an understanding of the recoil motion of D-particles, and thus of the 
associated Finsler-distortions of the neighboring space-time. For late eras, of particular relevance for us, the FRW 
Universe obeys a power- law expansion, this issue has been discussed in detail in [37|. In this case the deformation 
(|2.5p can be found by considering the geodesic motion of the D-particles. 

For completeness of our discussion let us review briefly the derivation of the induced metric deformations due to 
recoil in such a FRW universe background. We shall work with space-times of the form 

ds 2 = -dt 2 + a(tf(dX 1 ) 2 with scale factor a(t) = a t p , p < 1, (2.7) 

characterizing late eras of the Universe (radiation and matter dominance) of interest to us here. The pertinent geodesic 
equations in this case read: 

t + jrf 2 *- 1 ^*) 2 = 0, 

y + 2^{f)i = 0, (2.8) 

where the dot denotes differentiation with respect to the proper time r of the D-particle. 

For a single D-particle defect in the expanding Universe, for which the impact with a matter string state occurred 
at the time moment to, one has the initial conditions y l (to) = (for simplicity), and dy l /dt(to) = u l , with m the 
recoil velocity (|2.3p . In this case one easily finds that, for relatively (but not cosmologically) long times t ^> to after 
the event, the solution is: 

V l (t) = j^- (V- 2 ^ _ to ) + Op- 4 * ), t » t . (2.9) 

To leading order in t, therefore, the appropriate vertex operator describing the recoil of the Z?-particle in the FRW 
universe (|2.7[) is: 

V =f alT^-7T®(t-h)(ttl P -tot 2 Ad n X\ (2.10) 

where Q(t — to) is the Heaviside step function, expressing an instantaneous action (impulse) on the D-particle at 
t = t Q . 

To find the induced metric distortions due to the D-particle recoil (|2.10p . we first pass into the T-dual Neumann 
picture for the (spatial coordinate) X 1 u-model field (assuming T- duality to be an exact symmetry of string theory) 
of D-particle recoil and rewrite (|2.10[) as a world-sheet bulk operator [37| : 

V= [ g l3 y\t)d T X l = f d a {y i (t)d a X i ) = f ( Vi {t)d a td a X' + y l d 2 X t ) , (2.11) 

where the dot denotes derivative with respect to the target time t, a is a world-sheet index, and d 2 is the world-sheet 
Laplacian. This deformation is non-conformal, and one should dress it with a Liouville mode, which should then 
be identified with the target time. The discussion parallels the flat space-time case outlined above and will not be 
repeated here. Notice that it is the covariant vector yi which appears in the formula, and incorporates the metric gij , 
since y l = g^y 3 . 

To determine the background geometry, for the string motion, it suffices to use the classical motion of the string, 
described by the world-sheet equations: d 2 X % + T l ^ v d a X p 'd a X 1 ' = 0, where fi, v are space-time indices, a = 1,2 is 
a world-sheet index, d 2 is the Laplacian on the world sheet, and i is a target spatial index. The relevant Christoffel 
symbol in our FRW background case (|2.7p . is T 0i , and thus the operator (|2.1ip becomes: J s (jji — 2yi(t)T l oi) d a td a X l , 
from which we can read off an induced non-diagonal component for the space time metric 



2g 0i =y i -2y i (t)Tl i , 



(2.12) 



7 



In the FRW background, the covariant path yi{t) = gijH 3 {t) is described by: 

W(t) = ^ {< - to**) (2.13) 

(with ao the current value of the scale factor), which gives 

2g 0l = a 2 {t )u t . (2.14) 

This yields for the metric line element (after the identification of the Liouville mode with time): 

ds 2 = -dt 2 + u i a 2 {t )dtdX i + a 2 {t){dX i ) 2 1 for t > t . (2.15) 

As expected, this space-time has precisely the form corresponding to a Galilean-boosted frame (the rest frame of the 
D-particle), with the boost occurring suddenly at time t — tg. 

Such metrics can be generalized to describe the distortions of the space-time as a consequence of collisions of both 
open and closed strings with D-particle defects. Open strings are of particular importance to us, since they describe 
matter excitations of brane worlds, which can propagate in bulk space-times punctured by D-particles. As the brane 
world moves in the bulk, the D-particles cross it, and from the point of view of an observer on the brane they appear 
as space-time defects which "flash on and off". "Foamy" space-time features on the brane result as a consequence (cf. 
fig. [!}. When considering a collection of uniformly distributed D-particles, there are moments to, at which matter 
strings collide with defects. In (|2.15[) the time t is taken to be after the moment of impact but not much after on 
a cosmological scale. In a coarse-grained approximation, for non-clustering D-particle populations, one may assume 
that a(to) ~ a(t) to a good approximation. From now on it will be assumed that matter with momentum p and 
energy E induces the following space-time distortions (on our brane world, cf. fig. [TJ) as a result of its interaction 
with the D-foam: 

ds 2 ~ -dt 2 + Ui a 2 {t) (dtdX 1 + {dX 1 ) 2 ) . (2.16) 

Notice that in this formalism the time is already of Minkowskian signature, as a result of the Liouville dressing. 

The reader should notice that the recoil velocity of the point-like D-particles has been taken to have non-zero 
components only on the brane world, where the matter excitation resides. This is an approximation, which proves 
sufficient for the phenomenology of our model. Moreover, for the case of D-particles represented by wrapped-up 
D3 branes in type IIB strings [11|, the situation is complicated; the defects from such constructions entail additional 
degrees of freedom. For simplicity we shall ignore these complications, and concentrate on effectively point-like defects, 
whose recoil induces a back- reaction on space-time represented by a metric of the coarse-grained form (j2.16[) . 

In phase space, for a D3-brane world, the function m, (12.31) . involving a momentum transfer, Aki, can be modelled 
by a local operator using the following parametrization [38(: 

Ui = 5s 2AT = nki ' n ° sum * = 2 ' 3 ' ( 2 - 17 ) 

where the (dimensionful) variables fj,i = 1,2,3, appearing above, are related to the fraction of momentum that is 
transferred at a collision with a D-particle in each spatial direction i. These parameters are taken as gaussian normal 
random variables with a range — oo to +oo and defining moments 

(n) = 0, (2.18) 



(r irj )=0, iii^j (2.19) 

and 

a 2 = (r 2 )-(r. l ) 2 = (r 2 )^0. (2.20) 

A homogeneous foam situation will be assumed for cosmological eras after freeze-out of species, the situation of interest 
to us. It is, of course, understood that one may encounter inhomogeneities in more general models of D-foam, due 
to variations in the population of the D-particle defects at various epochs of the evolution of the Universe. However, 
such inhomogcncous contributions would be severely constrained by cosmological data at late eras of the Universe, 
and, in view of the weak effect of the foam, will be ignored in our analysis, to a first approximation. 

An isotropic foam situation would require r$ = r, for all i — 1,2,3. In such a case the variances (rf) = (Tq along 
all spatial directions. In what follows we keep our discussion more general, and treat r$ as differing along the various 
spatial directions. These differences are small, and could be possibly responsible for small-scale anisotropics in the 
large scale structure of the observable Universe, which, to a large degree, is isotropic and homogeneous. 



B. Supersymmetry breaking and implications for D particles 

An important issue [10( concerning our analysis pertains to effects of SUSY breaking in a phenomenologically 
realistic cosmological model of D-foam cosmology jj| E| ■ The important constituents in our construction are D- 
particles. The word particle in the latter is somewhat misleading. BPS D-branes in general ( and DO-branes or D 
particles in particular) are described within a conformal field theory as hypersurfaces on which open strings end. They 
can also be described as classical solitons (or defects) of an appropriate supergravity. For BPS D-particles (which 
break half the space-time SUSY) gravitational attractions are canceled by gauge repulsions (i.e. a no force condition 
holds). In this way, their gravitational effects in the expansion of the brane Universe, along the longitudinal (brane) 
directions, are suppressed. Typically the conformal field theory description is appropriate for weak string coupling 
(relevant for us) whereas the soliton representation is appropriate for strong coupling. However, in constructions of 
physical relevance, where many D-branes (through stacks, intersections and foam) are involved, a solitonic (i.e. defect) 
description is also possible at weak coupling. However SUSY is broken and the BPS condition will not hold. We will 
discuss the degree to which the no force condition still applies. 

Sen has shown that D-branes that break SUSY can be constructed in terms of branes and anti-branes (which can also 
be interpreted in terms of K theory) [39] . In fact the D- brane spectrum of some orbifolds of toroidal compactifications 
of type IIA/B superstring theory are compatible with the predictions of K theory. In general an orbifold of type IIA 
and IIB are generated by g\ = I n which denotes the reflection of n co-ordinates or g^ = 1n{— l)^ where Fl is the 
fermion number of left moving particles. When n is even the gi represent symmetries of type II theories. For n = 4a 
boundary state \D0) for a stable non-BPS D-particle has been analysed at a perturbative level in the string coupling, 
at tree-string level in ref. [40| and at one-string- loop level in ref. [4J| ■ The interaction of two D-particles separated by 
a distance r is determined by the amplitude 

(D0\V a \D0) 

where T> a is the closed string propagator 

_a_ f fz , Lo _ a 
° ^ii<i\ Z \ 2 

with a = ^ in the untwisted NS-NS sector and a = in the twisted R-R sector. The quantities Lq, Lq are the 
appropriate Virasoro operators. A similar calculation can be done for two D-particles in relative motion with a speed 
v. This analysis results in the absence of a force between two non-BPS D-particles for short distances up to order 
v 2 (the next order being 0(i> 4 )).We should emphasize that the no-force condition among the non BPS D-particles 
is found to occur at large scales r, as compared with the string length l s = \fa' \ where (target space) effective low- 
energy string action methods are applicable. This feature is exclusive of a critical-radius orbifold compactification, 
and points to the fact that specific compactifications may indeed preserve such conditions. However, the construction 
of [40] ignored higher-string loop effects. The incorporation of such effects results in general in the destruction of the 
no force condition |4lj . although under some circumstances it may be made valid up to one loop. Non-perturbative 
effects from string-loop resummation cannot be answered at present in this compactification framework. The above 
considerations will be used in section IIV A| when we discuss the contributions of the D-foam on the expansion of the 
Universe and its vacuum energy budget. 

There are other compactification schemes. Currently it is believed that phenomenology requires a spontaneously 
broken SUSY reducing to a N = 1 matter sector. This is hard to arrange in string theory model building. In models 
involving toroidal orbifold compactifications (as discussed above), attempts to obtain a realistic SUSY breaking 
scenario rely on non-perturbative effects such as gaugino condensation. The incorporation of D-branes altered the 
situation since they are sources of RR fluxes, permitting models with background fluxes. In our model we have 
a compact bulk dimension and fluxes cannot be turned on at will since they contribute to the vacuum energy. 
Orientifold planes are negative tension sources and are necessary to incorporate for consistency. For D-foam, we have 
however considered [42j a scenario of supersymmetry obstruction [43] rather than breaking, in the following sense: 
one may compactify the extra dimensions (relative to the three longitudinal directions) of the D8-brane world in the 
model of f o | 7 [Toll.f or the D7 brane worlds in the type IIB string model of using certain "magnetized" manifolds 

{e.g. tori) [44|. The resultant vacuum of the theory is still supersymmetric whereas the spectrum of excitations on the 
brane world (with three uncompactified longitudinal directions) is not. The bosonic and fermionic excitations couple 
differently to the "magnetic" fields in the extra compactified dimensions. Fermion-boson induced mass splittings are 
induced and can be made realistic [H| by appropriate choices of the "magnetic" field intensities -this is the analogue of 
the Zeeman effects in quantum mechanics. Thus the D-particles can be viewed as defects in a supersymmetric vacuum 
of a theory with a non-supersymmetric excitation spectrum; hence the no force condition would still be valid when 
there is no relative motion. Local fluctuations of the recoil velocities of the foam, as a result of the interaction with 
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stringy excitations, would break further the SUSY, as a result of the relative (slow) motion of the heavy D particle 
defects, but this phenomenon would not be sufficient to lead to significant gravitational effects of the D-particles 
on the brane worlds. Moreover, at late eras in the expansion of the Universe, the motion of the brane world in 
the bulk may be viewed as adiabatic and very slow, so again the relevant contribution to the breaking of the bulk 
supersymmetry would not be strong enough to create significant effects in the Brane Universe expansion along its 
longitudinal directions. 

Before dealing with the stochastic fluctuations of the cosmological space-time (|2.16p . we find it instructive to discuss 
aspects of "fuzzy" particle statistics in the current era. The evolution of the scale factor is ignored, i.e. the recoil of 
the D-particle is considered in a Minkowski background. This is done in the next section, and will help the reader to 
clarify the most important issues in our approach. 



III. PARTICLE STATISTICS AND THE D-FOAM MODEL 

At the macroscopic level there is some evidence [45j that the velocity distribution of self-gravitating collisionless 
particles follows non-extensive statistics [31] and is consistent with properties of dark matter halos in galaxies. Here we 
will show that, at a microscopic level, our D-foam model shares some features of non-extensive statistics for elementary 
particle dark matter candidates. The usual formulation of non-extensive statistics is based on the postulates of a 
generalized entropy S q and a (q) expectation value, (q) , where q is real. Explicitly, 

s q = 1^EL1$} (3.1) 

for a system with w microstates having a probability p t of being in the i-th microstate. Moreover we still require 
yi—i Pi = 1. As q — > 1 we recover the usual Shannon entropy. The q expectation value (0} q for an operator O whose 
expectation value in the i-th microstate is Oi, is defined as 

w 

(0) q = Y,O tP l- (3.2) 

i=l 

These postulates lead to modified particle distributions which, in the limit q —> 1, reduce to either Bose or Fermi 
particle distributions. Via maximum entropy arguments the generalized particle occupation numbers are given by 

( n J>9 = n i 72 TToZT" . . m i / To— i I c ( 3 - 3 ) 



[l + (g-l) / 3( £j .- /i )]V(<z-i)+£' 

where Ej is the energy of a single particle state labelled by j and £ = — 1, 1 for Bose-Einstein and Fermi-Dirac statistics 
respectively [3l|. The consequences of such a g-deformation of the usual statistics for "radiation" and "dust" en ergy 
densities in the early universe as well as for relic densities has been already discussed when ^ — 1 is very small [30( . 
Indeed this analysis used the available experimental bounds to deduce the constraint \q — 1| < 2.6 x 1CU 4 . In this 
approach q remains a parameter which needs to be determined outside the framework. 

In the analysis ensuing from the space-time foam model, the paradigm of D-particle recoil and induced metrics is 
allied to a conventional statistical mechanics (with an unconventional metric background). The parameter that we 
introduce will measure the "fuzziness" of space-time induced by the interaction of D-particles with stringy matter. The 
resulting randomly fluctuating space-time on a flat, non-expanding FRW background has the following form ( (|2.17[ l): 



/ -1 


rifci 






1 










1 


V r 3 k 3 









(3.4) 



The reader should recall that the (dimensionful) variables r^, i — 1,2,3 are taken as gaussian normal random variables 
with a range — oo to +oo and defining moments as given by (|2.18[) . (|2.19|) . (|2.20l) . Any potentially unphysical values 
of momenta introduced by this approximate model are rendered insignificant by the smallness of the variance of the 
Ti (suppressed as inverse Planck mass ). 



A. Energy-momentum dispersion relations 

For the statistical mechanical analysis we will need the energy-momentum dispersion relation in this unconventional 
metric. One way this can be found is by analyzing a massive spinless scalar particle moving in this background. In 
an arbitrary gravitational field v the massive Klein-Gordon equation for a field $ reads as: 
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[g^d^du - m 2 ]$ = 0. 
If we expand this for a general metric where all terms contribute, we have: 

g 00 (d ) 2 + 2g°%d i + Y / 9 U (d l ) 2 + 2g% tj d i d j - 



$ = 0, 



(3.5) 



(3.6) 



where Latin indices run over 1,2,3; in the interests of a streamlined notation, the Einstein summation convention is 
adopted when we have two identical indices but not when we have more than two. For plane- wave solutions 



to r satisfies 



or equivalently, 



&(x, t) = 4>(k, Lo r ) exp(i(—u! r t + k ■ a?)), 

-V'/.w. + £ .</"/••: + 2 (</"/,,/,-.),. + ,„-' = o, 

i 

g^kj ± sj (gWkj) 2 - g 00 (g^hkj + to 2 ) 



mi 
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This can be further simplified on noting the form of the contravariant metric: 
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(3.7) 
(3.8) 



(3.9) 



(3.10) 



where B = J2 3 j= i r j k ]- 

Choosing from (|3.9p the positive-energy solution ((— ) sign for the metric signature that we have chosen) and keeping 
terms up to rf, we obtain the dispersion relation: 



^^-E^ + fS^' (3 - n) 

i j 



where we dropped cross terms of the form r^j for i ^ j (these terms would in any case vanish in the end because of 
(HUD), i.e. used that: 



4 2 

r 
j j 



(3.12) 



j J j 

and E is simply a notation and stands for the standard Minkowski energy: 



E = 



(3.13) 



We remind the reader at this point that we use the standard Minkowski metric to transform between contravariant 
and covariant vectors, i.e. 



k l = ki, k° = uj r = —kg . 



(3.14) 



When the expansion of the Universe will be taken into account in subsequent sections, the full Friedmann-Robertson- 
Walker metric will be used. 



11 



B. Distribution functions 



From the standard version of statistical mechanics, in equilibrium, at a finite temperature T, we have for the particle 
distribution function (in units h = c = 1) 



1 



exp(/3(w r - n)) + £' 



(3.15) 



where £ = +1 applies to fermions and £ = — 1 applies to bosons, (3 = y and uj r is given in (|3.1ip . Writing down the 
Taylor expansion of (|3 . 1 5[) with respect to the small parameters rj and keeping up to terms r - yields: 



1 



exp ((3 (E-n)) + t 

exp (/?(£-//)) 
(exp 08 (E-n))+£) 

cxp (2/3 (E-m)) 
(exp 08 (£- M ))+£) 




(3.16) 



In our case, since r was introduced as a random variable, we will need to perform an ensemble average on (n) 
This averaging entails integration over the statistical parameters r< through integrals of the form: 



1 



dr., 



exp 




(3.17) 



where Oj is defined as: 



(3.18) 



In our case, the variance of the recoil velocity depends on the string coupling g s , which itself may fluctuate. Now 
we have: g~ 2 = e~ 2 ^ , where (0) is the vacuum expectation value of the dilaton field cj>. In a stochastic space-time, 
the D-foam Q, there are in general non-trivial fluctuations of the dilaton field, and, therefore, of the string coupling 
itself. Hence, our knowledge of such fluctuations is very limited since it depends on properties of string vacua. In our 
analysis we simply assume that aj are also normally distributed with \Pj) = (Jqj- At this point we should also note 
that <jQj are taken to be different along the various spatial directions, for the sake of generality. However the observed 
isotropy of the Universe in large scales imposes such potential anisotropies to be quite small, consistently with our 
ansatz of smallness of the <jj . 

On performing both averages on (n) r , the result, denoted by ((«)), on second order in the statistical parameters 
<jQj is given by: 



1 



exp(/3(£;- Ai ))+£ 

exp(/?(£ - fi)) 
(exp(/3( J B- Ai ))+0 2 



2/3 2 £ ty 



+ 



(exp (/3 (E-u))+ty 



(3.19) 



We will adopt the notation ((. . .)) to generically denote the above double average. 

One can see that (|3.19p has the same general form as the distribution functions appearing in (3(|, where Tsallis 
statistics was applied. Indeed, in that case, the corresponding expanded formula for the number density distribution 
function is: 



1 g-l (P (e - /j,)) 2 eH'-ri 



(3.20) 
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where q is the non-extensive parameter accounting for the non-additivity of Tsallis entropies. 

Comparing with our case here, we find that there are two major differences. The first is that, unlike in (30j . the 
quantity E appearing in the first term of (|3.19|) is not the actual energy of the particle, except in a Minkowski 
background (cf. (|3.13[) ). Therefore, in our case, by setting the stochastic effects to zero (which, in the Tsallis cos- 
mologies [30|| , corresponds to setting q = 1 in (|3.20p ). leads us back to the standard result by eliminating the last two 
terms in (|3. 19|1 : the metric f|3.4|) . of course, is also reduced to the Minkowski form, so that the actual energy gets the 
standard form (|3 . 13|) . The second difference is that the extra term in (13.20p depends only on the particle's energy, 
whereas, in our case, there is an explicit dependence on the particle's 3-momentum as well, owing to the induced 
Finsler-like metric distortion. 

The similarity of our approach with that of [30] , lies in the fact that the corrections to the number density distribution 
have broadly similar structure; the second term is a perturbatively small deviation from the standard case. In the 
Tsallis cosmologies, this is a consequence of the (assumed) smallness of the quantity q — 1, while in our case, this is 
guaranteed by the smallness of the statistical parameters aoj i describing the stochastic fluctuations of the D-foam. 

In the above sense, our D-particle foam may be viewed as providing a microscopic framework for modified particle 
statistics in the Early Universe, in a similar (but not identical) spirit to the Tsallis cosmological models of [30j. At this 
stage we would also like to recall the approach of [29], where a further non-extensivity in the framework of Tsallis, has 
been discussed, as being due to quantum fluctuations of the string coupling g s . Now g s — e'", where cj> is the dilaton 
field, and (...) denotes a vacuum expectation value in target space-time. By including all possible quantum space-time 
metric and topology fluctuations in a fully (yet unknown) quantum gravity framework, g s becomes stochastic. Such 
effects are ignored, as sub-leading, in our present approach. In principle they should, at least formally, be included in 
a truly non-perturbative formalism of the very Early Universe in the presence of space-time foam. 



C. Number and energy densities 



On applying the basic formulae: 



(2tt) 3 



((n))d 3 k (3.21) 



and 



P = 7^3 J (K»<f*. (3-22) 

with ((n)} given in (|3.19[) and ((nui r )) calculated in Appendix A, we can derive the number density n and energy 
density p for both fermions and bosons in typical eras for the evolution of the universe. The integrals in relations 
(I3.2ip and (|3.22p are with respect to fc J , which in the non-expanding case satisfies (|3.14[) . In our work here we restrict 
ourselves to non-degenerate species, so the chemical potential contributions will be ignored, since r»/i, with T the 
temperature. 

As our primary interest is the effect of space-time foam on heavy, non-relativistic Dark matter particles, we will 
start with this case. After some straightforward calculations, whose details are given in Appendix A, we obtain for 
the density of heavy non-relativistic particles, for which m 2 3> k 2 : 

n= JLj {{n))d?k 

= g ( ' e -(m-^ + J^S (mVip-VA a 2 ( -e"^"^ _ J^^i™-^ + J5 ^(m-^ 
\2tt(3) 57r 3 / 2 V r J °\2 16V2 27a/3 / 

9 (/T 2 + m 2 ) (^] V2 (te-^-M - JL^e-Hrn-riA t (3 23) 



127r 3 / 2 V ^ ' u V (3 J V 4^2 

where (Tq = Yli °f an d £ = +1 (£ = ~ 1) applies to fermions (bosons). 

As one can see from (|3.23[) , the first term on the right-hand side represents the standard result for the non-relativistic 
number density of fermions and bosons, whereas the other terms, proportional to CTg, describe the foam-induced 
corrections. It is interesting to note at this stage, that the D-particle foam seems to affect differently fermions from 
bosons, as manifested by the different £ dependence of the foam-correction terms in the distribution functions. This 
also holds for the relativistic case, which is included below for completeness, and the details of its derivation appears 
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in Appendix A. For the number density we get: 



n b = r(3)C(3) + ^T & 5 ^C(5) I -— r(6) + -r(7) ) (3.24) 



for bosons, and 



n f = ^2 |r(3)C(3) + ^T^ 2 ^r(6) [-2^(6) + 7,(5)] + y r(7) [677(7) - 6*7(6) + ^(5)} ) ( 3.25 ) 

for fcrmions, whereas the corresponding relativistic energy densities are found to be: 



^4 



Pb = |^r(4)C(4)T b 4 + ^a 2 T b 6 (^r(6)c(6) - ^r(7)C(6) + ir(8)c(6)) = ^ 3& r b 4 + s^o 2 ^, (3-26) 



for bosons and 



Pf = |^r(4)ry(4)r / 4 + 9 -La\T) (^r(6)„(6) - ^(7)^(7) + ^r(7)r?(6) + ^(8)77(8) - jjr(8)r?(7) + |r(8)^(6) 

7 V 2 ^4 , ^ 18941tt 6 50841 127i 8 \ (7\ tt 2 4 793.32 _ 2 ^ 6 

-3/?> + -^o?/ "777^ 5— C(7) + "T^r = U ^7Tff/ T / + 9f—^v T f , (3.27) 



8; 30^'"^"°^ V^5120 ST^'^^OQ-J [sj ^ yf ^~" ^ 

for fermions. 

These equilibrium distributions will be important later on when estimating the number of degrees of freedom of 
relativistic species. Changes to distribution functions induced by the foam can have cosmological implications, such 
as modifications to the freeze-out point of species and their (thermal) relic abundances, which could lead to observable 
effects, depending on the size of the effects and the associated string mass scale and coupling. In the next section we 
therefore proceed to discuss such cosmological aspects of the model. 



IV. COSMOLOGICAL ASPECTS OF THE MODEL 

The full rigorous cosmology of the D-foam is a difficult and rather complicated subject yet [l3| to be analyzed in 
a complete way. It depends on the underlying microscopic string theory model, which is not presented here [H, 
Below, we shall only concentrate on one aspect, which however is interesting for particle phenomenology, namely 
the effects of the non extensive statistics of the D-foam on thermal dark matter relics. However, before embarking 
on a detailed discussion of the effects of the interaction of the Dark matter with the D-foam on the thermal relic 
abundance of the former, it is instructive to discuss first briefly the effects of the D-foam fluid on the energy budget 
of the Universe per se. In this way the special role of the D-particle defects as supersymmetric vacuum defects, and 
not ordinary dark matter contributors, will hopefully be elucidated. We do so in the next subsection. 



A. D-particles, the energy Budget of the Universe and the Hubble Expansion Rate 

Our focus, for the interaction of D-particles with cold dark matter, has been on the induced stochastic Finsler-like 
metric fluctuations that are induced. Although D-particle is a generally useful term it is of course a DO-brane and 
represents strictly a space-time defect. Consequently it needs to be dealt with using the methods of string theory 
when considering its contribution to the energy of the system. We have noted that, if there is relative bulk motion, 
SUSY is broken and there are non-trivial forces among the D-particles as well as between the D-particles and the 
brane world or orientifolds The resulting non-zero contribution to the energy is proportional to v 2 for transverse 
relative motion of branes with different dimensionalities, and to v 4 for branes of the same dimensionality (there is no 
contribution to the energy of a p-brane-world from motion of other branes in directions parallel to its longitudinal 
directions). There is also a dependence on the relative distance of the various branes. In particular, the interaction 
of a single D-particle, that lies far away from the D3 brane (D8-compactified) world, and moves adiabatically with a 
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small velocity v± in a direction transverse to the brane, results in the following potential 2 Q 



87ra' 



^°o%8-+^ i -7^^»V«'. (4.1) 



On the other hand, a D-particle close to the D3-brane (compactified D8), at a distance r' -C Va', moving adiabatically 
in the perpendicular direction with a velocity v^ hort will induce the following potential to it: 



yshort ™'«° rt ) 2 U9) 



We now mention that the analysis @ of this short distance behaviour also reveals a characteristic minimum length 

Train ~ y/v±l S (4-3) 

which for jji < 1 is much smaller than the string length l s = vo', a feature common in theories involving D- 
branes (46| . The short range force is an attractive one and D-particles can make bound states with the D-brane through 
the exchange of light strings. This process in perturbation theory has divergences which lead also to instabilities. 
These are cured by a non-perturbative treatment of string theory, which reveals that the string coupling exp ((0)) at 
these short distances (r < ! 8 ) is given by [4(| 

exp (2 (0)) ~ gl (l + 9s ) (4.4) 
V {r/Va'yJ 

where g s is the string coupling on the brane for interactions between stringy matter. Clearly this coupling is enhanced 
and t associated mass of the bound D-particle - g^77^ is much reduced, compared to the mass (M s /g s ) of the 

D-particles attached on the brane world and co- moving with it in the bulk direction 3 . Hence, there will be two 
components of the D-particle contribution to the energy on the brane: one due to a D-particle density on the brane (a 
form of dust, in view of their no force conditions, discussed in section III Bp and the other due to the D-particle bound 
states which will contribute a negative interaction energy on the brane. Hence the presence of both contributions to 
the energy density allows a larger density of D-particles on the brane without overdosing the universe. 

Let us explain this point in more details. To this end, we first assume that the bulk space expands cosmologically, 
with a scale factor a(t), only along the directions that are parallel to the longitudinal dimensions of our D3 brane 
world (appropriate compactification down to three spatial dimensions is assumed from now on for the model of @, 
without further details; this by itself is an important topic for investigation that depends on the underlying models, 
but hopefully the main features of our analysis, regarding the contributions of the D-foam to the Hubble expansion, 
will not be affected significantly by such details). Let also rif = [3 be the three-space-dimensional density of 
D-particles on the D3-brane world (with Aft the corresponding number of (non-interacting) D-particles on the brane 
and a 3 (t) the proper volume of an expanding Universe for a co- moving observer) and n& = a /fy q the corresponding 
density of bound D-particles (with A/j, the corresponding number of bound D-particles, assumed non interacting among 
themselves). In view of their dust state, the D-particles bound on our brane world, and comoving with it in the bulk, 
will contribute a term in the energy density of the brane Universe of the form 

Pmass ~ ~ Tlf . (4.5) 

9s 

Since, as mentioned above, the bound state masses are much smaller, suppressed by powers of vj_ (for distances 
T ~ r m i n given by Eq. (|4.3|l ). due to the non-perturbative effects on the string coupling (|4.4[) . their mass- induced 
contribution to the ground state energy on the brane world will be subdominant, leaving only (|4.2I) as the dominant 
contribution from each one of them. Assuming, for simplicity, that the effects of the nearby bulk D-particles are 
effectively described by a model where these defects are mainly concentrated on a thin D3 brane parallel to our brane 



2 For brevity, in what follows we ignore potential contributions induced by compactification of the D8 brane worlds to D3 branes, stating 
only the expressions for the induced potential on the uncompactified brane world as a result of a stretched string between the latter and 
the D-particle - the compactitication does not affect our arguments on the negative energy contributions to the brane vacuum energy. 

3 That the mass of such D-particles is M s /g s can be determined by considering their recoil during scattering with fundamental strings 
having their ends attached on the brane, and moving in directions parallel to the brane's longitudinal ones |34|. |35t . 
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world, at bulk distances r m i n away, we obtain the dominant (negative) contributions to the vacuum energy density 
of the brane world due to the nearby bulk D-particles (with l s — 1/M S ): 

ptZly do = n b V^ds (r ~ r mm ) ~ -n b ^ JvZ ■ (4.6) 

Then the cancellation of the mass effect (|4.5I) of the D-particles on the brane world and those due to the short-distance 
potential U short USD gives 

12n f f*^ 

n b j= . (4.7) 



We will now use some estimates from a phenomenological study [47} of the order of magnitude of v± . In general, 
for a distribution of nearby D-particles in the bulk, this denotes the average velocity of the distribution. To simplify 
things, for the (qualitative) purposes of this section, we may make the reasonable assumption that the bulk D-particles 
are almost stationary in the bulk, and only our D-brane world moves slowly in the bulk dimension of the models of 
D-foam presented in Q . This makes v± the velocity of the brane motion, which has been bounded phenomenologically 
in [47} by the requirement that the models agree with the slow-roll conditions for inflation (defined appropriately for 
this particular class of brane cosmology models in ref. [J?}, where we refer the interested reader for details). The 
relevant phenomenological bound is 



and, on account of (14.71). this leads to 



v\ < 1.48 x lCT 5 ^- 1 (4.8) 



60 

n b ~ —rif . (4.9) 

gJ 

Another important estimate we need to make for our cosmological considerations in this paper concerns the density 
of defects, in both the bulk space and on the brane Universe. As we discuss briefly in section |V1 constraints on 
the density of D-foam are obtained by comparing the (subluminal) refractive effects of the D-particle foam vacuum 
in relation to time lags in the arrival of cosmic 7— rays depending linearly on energy (as reported by MAGIC [Hj 
and Fermi [4!| collaborations). In such studies it was found [l(| Eo| that, for a uniform D-foam density, at least for 
late epochs of the Universe's history, an order of magnitude rif ~ 1CP 3 — 1CP 2 was sufficient to make the D-foam 
explanation of the effect viable, in particular concerning the bounds imposed by the Fermi satellite [49j (for more 
discussion see section IVl and references therein). Since g s ~ O (1), we see that n b ~ 60 rif ~ 6x (10~ 2 — 10 _1 ) < 1 and 
we have the required (for consistency of our perturbative treatment) dilute limit of D-particle foam. The situation in 
a semi- realistic model, such as that of is complicated since there are stacks of D-branes, orientifolds and brane 
intersections (leading to a complicated interplay of the effects of U short f (14.2p ) and U long ( Q4.1[l )). However, we will 
consider the implications of the broad possibilities of D-particles on the brane and bound states. This argument has 
shown that nf ~ O (l0~ 3 — 10~ 2 ) does not need to overdose the Universe. 

The next issue we shall be concerned with is the contributions of D-foam to the Hubble expansion of the Universe, 
namely to the equation of state of the D-foam fluid. As already mentioned in the introduction section of this article, 
in view of their point-like structure in space, D-particles may be naively thought of as ordinary (massive, of mass 
M s /g s ) excitations of the vacuum, thus playing the role of dark matter [23|. However, this is not the case in the 
D-foam models [8|, as becomes evident from the afore-mentioned properties and will hopefully become clearer from 
the following arguments. If nD3-foam(i) = n f(t) denotes the number density of the D-particle effects on a D3-brane, 
at any given instance of cosmic time, t, then according to earlier arguments, we would have the following cosmic time 
evolution: 



37«D3-foam(£) + 3i/n D3 _ foam (t) = --r —7— (4.10) 

at V at 

where V ~ a 3 (t) is the proper higher-dimensional (bulk) volume, expanding only along the D3-brane world longitudinal 
dimensions; from the point of view of an observer on the three- brane Universe,^ Af is the effective rate of bulk D- 
particles crossing the brane as the latter moves in the bulk space (cf. fig. [T|). For eras of the Universe long after a 
catastrophic cosmic event, such as collision of two brane Universes in the arrangement of fig. [IJ the brane motion 
in the bulk is adiabatic to a good approximation. Secondly we assume that in such eras a steady state situation 
arises the effects of the source term in (|4.10[) compensate those of the cosmic dilution on the effective population of 
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D-particles in the expanding D3— brane Universe. Hence nD3-foam(i) , in (14. 10[) is approximately constant in time, 
i.e. 



In such a case, the D-foam would behave as a dark energy (cosmological constant-like) contribution to the energy 
budget rather than dark matter, thereby leading to different considerations than in the work of [2(3], where, as 
mentioned earlier, the D-particles are viewed as dark matter due to their localized nature in space. As the D3-brane 
moves in the bulk, such phenomena are obviously cosmological era dependent, given the likely inhomogeneities in the 
bulk population of D-particles. 

To summarise, the negative contributions to the Brane Universe's vacuum energy induced by averaging (14.21) over 
appropriate populations of bound states may then cancel out (or suppress significantly, on average) the positive energy 
contributions of the D-particles trapped on our world, due to their mass M s /g s . According to the above picture, the 
effective Hubble rate H(z) of the expansion of the D3 brane Universe, is not affected in the eras where (|4.11[) is valid 
by the presence of D-foam, thereby following the rule 



where z is the cosmic red-shift, f2, are various energy densities (in units of the critical density), and the . . . denote terms 
due to the brane Universes, such as dark radiation , which we do not discuss here explicitly. The quantity Qm contains 
dust contributions, including field-theoretic contributions to Dark matter (e.g. supersymmetric matter contribution 
coming from the effective (broken) supersymmetric field theories of the excitations on the D3 brane Universe) . On 
the other hand, JIde corresponds to dark energy/cosmological-constant-like contributions, with equation of state 
approaching w — > — 1 for redshifts z — > (today), which remain constant in time and dominate over matter for more 
than 1.9 billion years, according to observations. 

Although a detailed and complete microscopic bulk model, with realistic particle phenomenology and Cosmology, 
is far from being available in this framework at present, nevertheless, we consider the above picture as plausible, 
and this will be our main assumption in what follows. For a more detailed, but definitely incomplete, microscopic 
picture, where one takes into account also the bulk distributions of D-particles lying far away from the brane world 
(cf. Eq. (|4.ip ). we refer the reader to the second paper in [10|j . The possibility of approximately constant D-foam 
contributions to the vacuum energy (|4.11[) was also raised there. 

Based on the above assumptions and plausibility arguments, we shall assume in what follows that the dominant 
effects of the D-foam on the energy budget of the Cosmos concern primarily the (ordinary) Dark Matter content of 
the brane Universe. These are in the form of the distortion of the neighbouring space-time during the propagation 
of (electrically neutral, weakly interacting) matter in regions of space-time with populations of D-particle defects; 
this leads to an effective background space-time of the form (|2.16p . felt by the Dark Matter particles. This space- 
time fluctuates in the framework of stochastic D-foam according to (|2.17p ). (|2.18p . (|2.19p , (|2.20p . Such space-time 
distortions depend on the momentum transfer between the defect and the matter string (i.e. akin to Finsler metrics), 
and as such affect the relic abundance from the Boltzmann equation for the evolution of Dark Matter. We now 
proceed to discuss this effect in some detail, in the context of thermal relics that we restrict our attention upon in 
this article. 



We now come to the announced detailed discussion on the effects of the D-foam on thermal Dark Matter abundances, 
as a non-trivial consequence of the distortion of space-time (|2.16p on account of the non-trivial interaction of the 
pertinent matter string excitations, representing Dark Matter relic particles, with the D-particle defects. To this 
end, we shall first need to calculate the effective number of degrees of freedom participatin g in the thermal history 
of the D-foam universe. The reader is invited to compare our findings with the approach of [30], where conventional 
cosmology but with Tsallis statistics [3l| has been considered. If one assumes that the dominant energy contribution 
comes from the relativistic species (e.g. the radiation dominated era (r.d.e.)), then the total energy density can be 
written in the very compact form: 



-T7«D3-foam(i) — . 



(4.11) 




(4.12) 



B. Effective number of degrees of freedom 



PT = ^9effT 



(4.13) 



17 



where g e g is called the effective number of degrees of freedom and in the standard approaches is given by: 

4 _ /rr 4\ 



9eff- 



where Tij, (T^t) denotes the temperature of the thermal distribution satisfied by bosonic (fermionic) species i; in 
principle, ^ T, where T is the temperature of the photon gas. The quantity (I4.14[) is very important in cosmology 
because it allows us to count the species that are relativistic at a given era. Obviously it is a time (temperature) 
dependent quantity since, as the Universe expands (and freezes), less and less species can be effectively massless, i.e. 
highly relativistic. For instance, at very early eras, when T > 175 GeV, all the present species were relativistic yielding 
a very big effective number of degrees of freedom: g e g ps 106.8. After the e~ annihilation, at T > 500 keV, only the 
three light neutrinos (restricted to one helicity state) and their anti-neutrinos and the photon remain massless, giving: 
g e ff~ 3.37 (for more details on this see [HI)- 

Within the context of our model, the effective number of degrees of freedom is expected to be modified. To obtain 
this correction, let us substitute in the left hand side of f|4. 1 3[) the relativistic energy density for fermions and bosons 
that we derived in the previous section (expressions (|3.26[) . (|3.27|) ) . Then we have: 

30_ 2 /2^ 4 ^ - fT hb \ 4 793.92 ^ , fT u \ 4 \ 

Now the cosmological evolution equations, for instance in a r.d.e, will still be given by the standard relations: 

H = lM(g' eff )i— (4.16) 

" rn.pi 



and 



t = 0.30(9^)-^, (4.17) 



but with the effective number of degrees of freedom having the corrected form (|4.15p . that carries the D-foam effects. 
Hence, it can be seen that all the additional effects stemming from the D-particles recoil, can be entirely subsumed 
in a modification of the effective number of degrees of freedom as presented by (|4.15[) . 



C. Boltzmann equation 

In order to account for the expansion of the universe by the factor a(t) we should rewrite the distorted metric (|2.16[) 
in the following form: 




(4.18) 



Hence the relations between covariant and contravariant momenta now become: 

fc = -k°, h = a 2 (t)k\ i = 1,2, 3, (4.19) 

where in our notation the scale factor a (t) is dimensionless and the recoil effect has not been taken into account at 
this point. We also have the general relations: 

= m—,n = 0,1,2,3, (4.20) 

(XT 

where r is an affine parameter (the proper time in our case) and m is the mass of the particles under investigation. 

We should point out this perturbation is different in spirit (and mathematically) from metric perturbations in cos- 
mological perturbation theory. Our considerations of quantum gravity phenomenology suggest that an effective theory 
cannot be just based on a 4-d pseudo-Riemmanian manifold M but rather the tangent bundle TM. Diffeomorphisms 
would have to be appropriate to TM and hence issues of our metric being diffeomorphically equivalent to those used 
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in cosmological perturbation theory cannot materialize. To illustrate this, in the conformal Newtonian (longitudinal) 
gauge, a line element in cosmological perturbation theory can be written as 

ds 2 = a 2 (r)[-(l + 2*)dr 2 + (1 - 2<Z>)5 ij dx i dx j ] (4.21) 

where 'J and $ are functions on M. It is interesting to see how similar we can make the metric in (|4.18p to that in 

To this end, the reader should recall that the metric (|4.18|) has been derived by a Liouville-dressing procedure in 
the world-sheet conformal field theory describing the propagation of a string excitation in a background space time, 
as reviewed briefly in section III Al It can be diagonalized by means of an appropriate T duality transformation on 
the target-space coordinates/cr-model fields, as follows [50]: Ignoring for the moment the cosmological expansion, i.e. 
setting et(i)=constant, we observe that the vertex operator (12. ip for the description of the recoil of a D-particle can 
be rewritten, using the two dimensional version of Stokes theorem, to give [50] 

£ d 2 zu l e a pdh[Q e {t) + tS £ {t)]d a X l , (4.22) 

where a,/3 — 1,2 are two-dimensional world-sheet indices. For times long after the impact of the open-string state 
with the D-particle, t > 0, this can be interpreted as an open string propagating in a target-space antisymmetric 
background B^ v with 

Boi ~ ui and Bij = . 

On making a canonical T duality transformation of the co-ordinates, the presence of the B field leads to the following 
mixed boundary condition on the boundary <9E of the world-sheet E: 

where E is homeomorphic to a disc. The induced open-string effective target space-time metric is obtained from 
calculating the world sheet propagator (X ,1 (z,'z)X^ (0, 0)). The resulting metric is [5(| 

9nv = (1-^ 2 )V" M,^ = 0, 1 

9n,v = i)iiv, |ti> v = all other values , (4.23) 

where we have considered a frame in which the matter particle has motion only in the spatial X 1 direction. Also, it 
should be noted, the metric implies the existence of a critical recoil velocity (viz. the speed of light in vacuo), which 
stems from Lorentz invariance of the underlying string theory [50l | . 

On embedding such diagonal metrics in Cosmology, one should multiply the spatial parts by the appropriate 
Robertson- Walker scale factor a(t). However, the reader should recall that the dependence of the metric elements 
(|4.23|) on the recoil velocity implies in turn a momentum-transfer dependence of the metric. Then, compared to (|4.21[) . 
although in diagonal form, the space part of the cosmological extension of the metric (|4.23p has a velocity dependent 
piece, quite different from metric perturbations used in cosmology, that depend only on coordinates. Also, the time 
part has a similar form. Consequently, the Finslerian nature [2j of our metrics differentiate them significantly from 
usual metric perturbations. 

Because of the above differences, it is not possible to put the back-reacted metric (|4.18p into the longitudinal 
gauge of scalar metric perturbations by an ordinary (i.e. only a space-time coordinate dependent) diffeomorphism. 
Consequently, the metric due to the clumping of ordinary particles cannot be put into the form of the metric due 
to back-reaction from D-foam. As we have argued previously, the role of the D-particle defects is quite different 
from ordinary dark matter. As discussed above, in view of their argued no-force condition, which may be maintained 
in our supersymmetry obstruction scenarios characterizing models of D-foam, these defects constitute a dark fluid, 
contributing to the dark energy of the vacuum [10J. The fluid affects the propagation of ordinary particle excitations on 
the brane world through the above-mentioned Finsler-type metric distortions (|4.18p , and via this the relic abundance 
of ordinary dark matter particles. This will be discussed in detail below. 

The momentum dependence of the (Finsler) metric implies a modification of the geodesic equation for along 
with the expected modification of the Christoffel symbols. Details are provided in Appendix B. 

However the distribution function of a particle species / is specified by the microscopic properties of the system 
and, in this sense, it would be convenient to define a local, in space-time, momentum in an expanding universe (5l) : 



k l = a(t)k l , i = 1,2,3 



(4.24) 
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In terms of these re-scaled momenta, the energy-momentum dispersion relation for a (dark matter) particle of mass to 
in the spatially flat Robertson- Walker space-time background we are working on, assumes an effectively "Minkowski- 



space-time" form, k^k 1 'g^ v = —uj 2 + a 2 {t)k l k^ 8i 



-lo 2 + k^k J 6- 



-m 2 and hence it is essential to define the 



phase-space densities as functions of the coordinates {x l , t) and the local momenta k , f(x l , t, k ). This is particularly 
important in keeping the correct scaling properties of the Dark Matter density with the scale factor, as we shall see 
below. In an isotropic Robertson Walker background, the momenta are assumed to be on-shell, so that the phase-space 

density depends only on the amplitude |fc|, and thus the energy oj, of the Dark Matter particle. However, in our case, 
the small anisotropics that characterize the foam fluctuations should be taken into account, and this is performed by 
assuming a dependence on the individual components fc, as well. Then, the Boltzmann equation will be written in 

terms of the one-particle distribution function / ^r^fej and explicit dependence on the energy (jt° = k°^j can be 

omitted due to the on shell condition: 



g^k^k" = -to 2 . 

For such a distribution function, Liouville operator takes the form: 



L[f] = fc' 



df 



E 



df dk 



g x fi Z-j gj^t d T 

and after applying the isotropy condition: = for i = 1, 2, 3, it simplifies to: 



(4.25) 



(4.26) 



Noting that: 



t r n TP df i df dk 



dk , . dk . . . dt , 4 f . dk . , . k , ,- 
— = a(t)— + d(t) — k l = a(t)— + a(t) — k 
ut dr dr dr to 



and using the notation H — ^, we get the final expression: 



L[f] = k —+m - 



df 



dr 2 



dk 



dk 



(4.27) 



(4.28) 



(4.29) 



Substituting (|6.36l) in the above, with k° being everywhere replaced by k , and dividing throughout by k , we 
obtain: 



L[f] _ d£ 
dt 



H yV^l- 2Ha 2 (t)k° Y rj? ^- + 8ff a 4 (t) V r? (*T — 



df 



df 



k° 



dk 



dk 



dk 



±Ha 2 (t)y(VYyr-k^ + AHa\t) (VY^r^^ - ±a\t)H^ (fcT£rffc^ (4.30) 



dk 



df 



dk 



df 



3 % » 3 

Applying at this point the mass shell condition (I4.25[) . yields the energy dispersion relation: 

2 



dk 



k° = a 2 (t)J2(k 



r, 



to 



1 



At) (Ei(fc*) ! 



TO^ 



(4.31) 



Further analysis involves the approximation of heavy dark matter, which we consider in this work, since the phe- 
nomenologically dominant species are the heavy ones that in general are expected to cause the most significant 



-H 2 



distortions in the space time background, to 2 ^> J^. k , in which case 



1° 



in 



2to 2 



a 2 {t) 



rn 



(4.32) 
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or expanding also up to second order and dropping cross terms, 



1 1 



(4.33) 



With these in hand, (|4.30j) gives: 



(4.34) 



Now we can average over the ensembles for the random variables r^. It is not trivial why this can be done at the 
equation level. The reason is because we assume that the time scale of D-particle scatterings is much shorter than 
H . Remember also that terms r^j that would have zero contribution for i ^ j, have already been ignored, in 
previous steps of our analysis. Then, the Boltzmann equation (with the binary collision term taken also into account) 
is written as: 

taA'l^ifS - Wa'(t) •£ (WfY: »l$~,= ^1, (4-35) 

^ dk j 2 c^A: /c 

where the fourth and sixth terms can be neglected as very small compared to the other terms, always under the 
assumption of superheavy dark matter: ( — ) < 1. Since the number density n (t) is defined as: 



»(*) = _* / d d kf[t,k'), (4.36) 
(2tt) 



^ 2 3 

with d 3 k = dk dk dk , the following simple identities can be derived: 

->— d „ { —i\ dn 
d 6 k- 



f(t,k l ) 



(2tt) 3 J dt \' / dt' 



dl 



~3n, 



7-^73 ^2 a 0i f d 3 kk l ^L = - (<j 2 01 + a 2 02 + a 2 3 ) n 

(27T) . J dk 

We now define an average temperature T through (see [52]): 

^3 J d 3 k (fc') 2 / = Tmn (4.38) 



and, as a result, the last identity in (|4.37p becomes: 

9 2 f .37: /VV °J _ IJ2. 1 _2 ,2 



(2^) J Y ^ W dk 
Therefore integrating (|4.35p with respect to d 3 k yields 0]: 



3 E 4 fd 3 k(k l Y^l=- « + + o* 3 ) 3Tmn. (4.39) 
•/ K ' dk 



_ + 3ff „ = r(t ) „ + ^ s y^ I (4.40) 
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where we have managed to incorporate all the foam-induced correction terms in a time-dependent source term given 
by: 

r(t) = 2Ha 4 (t)m (cr^ + ct§ 2 + o% 3 ) [9T + 2m] . (4.41) 

The reader is reminded at this point that the above result has been derived under the approximation of superheavy 
dark matter species \- 

We will elaborate a little on the contributions of the foam to the interaction terms C [/] in the Boltzmann equation 
(|4.40p . We assume for simplicity that the dominant scattering process affecting the abundance of x is 1 + 2 < — > 3 + 4, 
where 1 and 2 are \ particles and 3 and 4 denote light (Standard Model) particles. The latter are assumed to be 
strongly coupled to the cosmic plasma and hence in equilibrium. Following standard treatments we have 

i r 3 c[f] r (Pk 1 r d 3 k 2 r d 3 k 3 r d 3 k 4 

Wf J "17 = J (2^) 3 2 Wl J (2^) 3 2w 2 J (27r) 3 2w3 J (2tt) 3 2w 4 

x (2tt) 4 (5 3 (k 1 + ~t - k 3 - fc 4 ) 5 ( Ul +lj 2 -uj 3 - w 4 ) \M\ 2 

x [hh (1 ± h) (1 ± h) - hh (1 ± h) (1 ± /*)] , (4.42) 

where (±) denotes Bose enhancement (+) and Pauli blocking (— ); fa is the distribution function associated with 
species v. Ai is the scattering amplitude; and u>i is the foam dressed energy for particle i. For T < uj — /i the quantum 
statistics of particles can be replaced by the Maxwell-Boltzmann distribution, i.e. fa — > exp (3 (ni — u>i) and clearly 

the effects of foam in fa will only occur through the modified u>i {= k°^j . D-foam effects will also appear in M.. For 

neutralinos (which are spin-i Majorana fermions in MSSM and are linear combinations of the bino and wino) the 
tree level amplitude involves a weak vector boson propagator. Vector bosons have masses close to 100 GeV, i.e. of the 
order of WIMPs and so (for simplicity working in the unitary gauge) the foam effects at the tree level are governed 
by the vector boson propagator; this leads to a suppression of the effect of the foam relative to that in the source 
by the square of the WIMP mass. From the form of the Boltzmann factor /, the thermal averaging has a similarly 
suppressed foam contribution. 

Finally, it is important to note that equation (|4.40p yields the standard scaling for ordinary massive matter in the 
absence of D-particles (term 3Hn) , because of the use of the local momenta formalism [5lJ . From (|4.40[) we can also 
see that stochastically fluctuating Finsler metrics, induced by the D-particles foam model, are equivalent to particle- 
production source terms in the Boltzmann equations. As a consequence of these extra terms, relic abundances of 
dark matter candidates are expected to be modified, more or less significantly, according to the free parameters of our 
model. To understand this modification, the reader is referred to |25| . where the generic effect of source terms on the 
Boltzmann equation is investigated, but where the source is generated by off-shell / dilaton terms. Another difference 
from [25| is that, in our case, the foam-induced source effects come with a positive sign, therefore are expected to 
increase today's dark matter relic abundances, when compared to the standard ones and as a result, leave even less 
room for supersymmetry at colliders |25|. However, the corrections will in general be quite small since they will always 
scale as , which are perturbative free parameters in our analysis depending on the string scale, as will be explained 
at the end of the article. 



D. Solutions to the Boltzmann Equation and Foam-Modified Thermal Dark Matter Relic Abundances 

In (I4.40p . from standard arguments, the collision term has the form: — (av) (n 2 — n^ q ) [28]. Then the Boltzmann 
equation is re- written as: 

h + 3Hn = t(t)n - (av) (n 2 - n 2 cq ) , (4.43) 

with n cq being the (thermal) equilibrium distribution, (crv) the thermally averaged annihilation cross section summed 
over all contributing channels and v the Mueller velocity. Our source, from now on, will be denoted by T(t), not to 
be confused with the particles interaction rate T(t). 

An analytical solution of (|4.43p is difficult to be given generally since the collision part is era dependent. In standard 
approaches [28], the total annihilation cross section is assumed to depend on temperature through: 

(av) = a(x) = cox~ n , (4.44) 

where x = ?S , Cq is an arbitrary constant (discussed in more detail later) and n is related to the channel of annihilation: 
n = for s-wave and n = 1 for p-wave annihilators. Since these are usually the dominant types of annihilation, we 
will restrict our attention on these two possible values for the integer n. 
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It proves convenient to use instead of the particles' number n(x), the parameter number per entropy density 
Y(x) = liipl. Under this re-parametrization and along with using the conservation of the entropy per comoving 
volume (sa 3 = const =>■ = —3Hs), (|4.43p is written in the form: 

dY 

— =T{t)Y-{av)s(Y 2 -Y 2 q ). (4.45) 

Also, since the interaction term depends explicitly on temperature rather than time, it is useful to take x as the 
independent evolution variable. Applying (|4.17[) then, yields: 



y_^0± {Y *_ Yelh (4 . 46) 



where we have set w = ^or' 1 m3 > with h denoting the entropy degrees of freedom and H m = Hx 2 = 1.67 g'^g 2 jj^, hi 
a r.d.e, with H being the Hubble rate. 

Before the freeze-out (x < Xf), particles are in thermal equilibrium, therefore 

Y{x) « Y cq {x) (4.47) 

with the equilibrium density being: 

Y cq (x) = n eq (x)/s = u ~ 1 a; 3 n oq (x) (4.48) 

and n eq (x) is given by (|3.23p . 

After the freeze-out and at relatively late times (x 3> Xf), Y tracks Y cq very poorly, therefore assuming that 
Y ^ Ycq, one can safely neglect Y eq from (|4.46p obtaining: 

d 1 ! x F(x) 1 = (cry) up 1 
Hm x Y X 2 

where we have also divided both sides by l/Y 2 and have restricted the discussion on species \ t° superheavy dark 
matter candidates of mass m x . Now (|4.49|) is in the form of a Bernoulli equation (y' + P(x)y = Q(x)y n , for n = 0) 
and as such can be exactly solved. The general solution to (|4.49D is then: 

with a and b two arbitrary constants and s just an integration dummy variable, not to be confused with the entropy 
density. Choosing a = Xf, keeping b general and applying (|4.50[) on the freeze-out point Xf as well (this is only 
approximately correct since the solution (j4.50p is formally valid for x > xf), we can obtain today's solution: 

- ex P (- r x -M d \ i = exp (_ r *jM dx ) S r exp ( r ?m dx ) ^ ds \ , (4 . 51) 

Y{x Q ) l \ 4 H mx )Y(x f ) P \ J Xf H mx )\J Xf P \J Xf H mx J H mx s 2 f> { ' 

where obviously Xq denotes today. Expanding the exponentials with respect to the small parameters erf that are 
hidden in the source term T(x) and keeping only corrections up to order of, that is terms linear in T, one obtains: 



Y(xo) Y(x f ) 

' " #H + r r i^) + r ( /' M ^ 

f n m x J * \ x f) Jx f n m x ° \ J Xf £1 rn x J J Xf \Jxf T1 rn x J J^m x o 

The function cr(s) should in general also depend on the statistical parameters oi of our model through cq. However, 
at this point this is not taken into account since formally the only effect it would have is to add extra undetermined 
constant terms in our analysis. An important point to be remarked here is that the freeze-out point Xf, is also affected 
by the source term (|4.41[) . To get the shift to Xf we apply the freeze-out criterion [28| : 

Y(x f )-Y oq (x f )^dY cq (x f ), (4.53) 
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where d is a phenomenological constant. Using (|4.50[) and (|4.48[) . we get: 

xT(x) d \ ^fH ds = (d+i)-^- 3 ^- 1 ^/), (4.54) 

with n e q(xf) given by (|3.23p . 

Denoting the left hand side of ()4.54[) by I, a very rough estimation of the freeze-out temperature comes by assuming 
that: I w /(°) where by /( 0) we mean the source-free part. 

Then, choosing b = x in , where a^„ stands for some initial time, e.g. after just after the exit from inflation and 
expanding the exponential with respect to the source T(t), we get: 
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since Xi n <C X/ and _B is already of order <r 2 (see below) . We also get 



(4.55) 




I (0) = -At \-iLt--L), (4-56) 



n + l \ T (°)" +1 
\x f 

m 2 

where A = and 5 = 18-^CTq. Applying our assumption then and always to linear order in B, yields: 



(0)\ "+ 1N 



or equivalently 




n+l 

2 =.2 



/n (Jq. (4.58) 



We can see from (|4.58j) that the freeze-out point is positively shifted due to the foam meaning that dark matter species 
are expected to have decoupled later in the presence of the stochastic background. However what is important is that 
the correction to Xf, similarly to [25J, only scales linearly to the source term, therefore as of. And since in any case 
these parameters are very small, also the shift on the freeze-out temperature will be small, therefore one can retain 
the assumptions/approximations of Standard Cosmology [25|,[28j]. In this sense, all terms involving 1/Y(xf) can be 

neglected from equation (|4.52[) . For convenience, we additionally can set f*° ?jr-ds = J, as in standard approaches 

and = J( s ): where by xj. We note at this point that from now by Xf, wherever it appears, we will refer to the 
modified freeze-out point. Then, equation (|4.52[) can be written under the compact form: 

Yixo)- 1 = (J) m x M pl g'-Jfh'J (l - £ ^dx + i £ J (a) ( jf <fe) , (4.59) 

where g' e g * and h! represent the effective number of degrees of freedom and the entropy degrees of freedom at the 

freeze-out point Xf, respectively. Obviously, the term (j^) 1 ^ 2 m x M p i g~g^h! J corresponds to the standard result [28j 
whereas the other terms in the prefactor are our foam-induced corrections. However, one should bear in mind that 
since both the freeze-out point and the effective number of degrees of freedom carry foam corrections, what we refer 
to as standard result is only approximately the standard one. We also note, that solving exactly the differential 
equation ()4. 59)) gave us an extra correction term (third term in the prefactor in (I4.59[) ) when compared to [25j, where 
similarly a source term in the Boltzmann equation was considered. 
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The energy density of the species \ wm then be: 

r xo xf(x) 



P = Y(x )s 



or equivalently to [25]: 



1 



-dx 



J(s) 



3 t* g' 



7 VeffJ 
M p l J 



(4.60) 



(n x h%) 



no source 



9effJ 



x ° xT(x) , 1 

~H dx -- T 



J(s) 



xf(x) \ 



ds 



(4.61) 



It is important to note that the formalism presented so far is true for an any source term present in the Boltzmann 
equation, only under the restriction that it is perturbatively small. Substituting in then our specific source term 
(|4.4ip and using relation (|4.44j) for a general n, the integrals of the first factor on the right hand side of ()4.61j) can be 
explicitly calculated yielding: 



x ° xT(x) , 1 
~H dX ~ 1 



J(s) 



xT(x) 



1 + x A m 2 x al 



t 

6xq Qx j ~\~ Xq x j —j~ 



dx }ds> = 



„2+n 



ds + 6 



-1+n 



ds 



with n — 0, 1. 

Upon noticing that in all realistic cosmologies [53j <C 1, 



-2+n 



ds 



(4.62) 
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n + l\ x n t +1 



„n+l 



n + 1 x n t +1 



and therefore for both values of n, the following approximations can be made: 
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(64 + 4) 



C X 



,2+ri 



J 



cpx 

J 



1 



Xf 



s -2+n 
1 

„— 1+n 



ds 



ds 



ds 



6x~ n 



n+1 



EL 

Xq 
n + 1 



ra+1 



■Li 

Xq 

n+l 



6Xn 



« 1, 



4+n 



XQ J 



El) 

XoJ 



5+n 



«1, 



< 1. 



Hence, f|4.61[) reduces to: 



Cl x h 



{tt x hl 



X' L J no source 



9 effJ 
9effJ 

9effJ 
9effJ 



■'■ 



-n+l 



„2+n 



(l + m^o^l + exo- 1 )} 



(4.63) 



(4.64) 



(4.65) 



since terms of the form (^r^j are of subdominant contribution. It is important to realize that no era assumption was 

made when calculating the integral f*° j^^- dx in (|4.6ip although the form of the source term as written in (|4.4ip 

assumes a r.d.e. In fact, one can easily check that the integrand does not require any specific form of the Hubble rate 
H and thus is era independent. 

However, totally in our calculation, it has been tacitly assumed that temperature scales as T = throughout 
all the evolution of the Universe. This assumption is only approximately true, since annihilating particles at different 
stages of the evolution of the universe deposited energy in the universe, slowing down its cooling. In order to get the 
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corrected cooling law one should assume that the entropy, dominated by the relativistic species contribution, remains 
constant (this is also an approximation since the presence of the source term may lead to entropy production [25]). 

1/3 

Then we get: a(t) f /ffit) ) Ss t l i where g e g is an evolving with time (temperature) function g e ff(t) [28[ and is such 

that today g e ff,o — 3.36 (only photons and the light neutrinos contribute), while at freeze-out g e ff(tf) — 106.75 within 
the framework of the Standard Model. Therefore in principle, one should calculate the integral in (|4.61[) numerically 
using this evolution of g e a(t(T)). For our qualitative purposes in this work we have ignored such corrections which, 
at any rate, could not affect significantly the order of magnitude of foam contributions to DM relic abundances. 

The last ingredient in our calculation involves the modification of the effective number of degrees of freedom, which 
has been presented in section IIVBI Applying (14. 15)) and counting carefully the relativistic species at an early era of 
the universe, when basically all species of the Standard Model were relativistic (the reader should recall at this stage 
that a typical freeze-out temperature is of the order of some GeV for typical dark matter candidates with masses in 
the range m w 1 GeV — 10 4 GeV), we find: 

g' eff = 106.75 + 22138 afT 2 (4.66) 

(we used that: g f = Squarks+flieptons = 6x(2 x 2 x 3)+3x(2 x 2)+3x2 = 90 and g b = .9 g iuons+.9EW+.9 P hoton+gHiggs = 
8x2 + 3x3 + 2 + 1 = 28). 

Therefore the correction to the effective number of degrees of freedom will be roughly going as: 



,9 eg 

and the modification to dark matter relic abundances is finally found to be 



— ) » 1 + 207.38 a 2 Tf (4.67) 



(n x h 2 ) 



no source 



{ 1 + m 2 , a\ (6xq 1 + l) }{ 1 + 103.69 a\m\ xj 2 } . (4.68) 



The dominant contributions come from the terms (J^m 2 , and 103.69 a^m^x ^ 2 , since the freeze-out for, e.g. neutralino, 
dark matter candidates typically occurs at Xf ~ 20 [53J (since any foam- induced freeze-out shift has been found to be 
negligible, we have retained the Standard Cosmology estimate). The other terms can be safely neglected from (|4.68p 
since xj <C x$ in all phenomenologically realistic cosmological models. 

Then, writing our statistical parameters <7oi as: 

a Ql = ^-Ai, (4.69) 



where g s is the string coupling, M s the string scale and Aj are dimensionless variances ^ \fTJrj)) = A,, i = 1,2,3 
that can be naturally up to 0(1), the dominant contributions of the foam effects on DM relic abundances is of order: 

— ^1 1 + 1.259^-| £ A,?. ( 4. 7 0) 

V X 0/ no source JVJ s ^ 

Then, we may use the standard WIMP estimate (f2 x ) no SOU rco ~ p f M P (^A^) -1 , where p c is the critical density of the 
Universe, To is the CMB temperature today, and the total annihilation cross section can be written on dimensional 
grounds as: 

a A v ~ k-^f^ , (4.71) 

for s-wave annihilators (cf. (14.441) ). that will restrict our qualitative considerations here. The weak interaction coupling 
constant g woa k ~ 0.65, and the fudge factor k [53j, which in general parametrizes deviations from the above estimate, 
is assumed in the range k G (0.5 — 2). As well known, the WIMP miracle then implies that, for a particle which is 
assumed to make up 100% of DM, as is our case here, one must have masses in the range m x ~ 100 — 10 3 GeV. In 
view of (|4.70p . then, we obtain (for s-wave annihilators): 
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We shall make use of such estimates in our phenomenology analysis, in the next section. 

It is understood that in estimating J7 no sourco one may use the pertinent thermal abundances in the context of specific 
particle physics models of interesting dark matter candidates, such as supersymmetric or string-inspired models, with 
hidden sectors [53]. The Finsler foam corrections are universal in this respect, since, as we have seen above, the only 
assumption is that electrically neutral candidates interact predominantly with the foam defects. The only difference 
in using models other than the standard WIMP ones lies on the fact that in such cases the weak interaction coupling 
constant g wca k should be replaced by the appropriate coupling gx , for instance describing the properties of the hidden 
sector of a string-inspired model [53]. The restriction is that the combination m x /g x , which enters the inverse total 
annihilation cross section of the dark matter candidate X, is fixed by the requirement of producing thermal relic in 
the ball park of WMAP observations [13. [l^. 

The presence of the stochastic D-particle foam, results in an increase of the thermal DM relic abundance, compared 
with the corresponding result in Standard Cosmology, in the absence of the foam. In our approach, the modifications 
are mainly due to the presence of a source- like term (|4.4ip in the pertinent Boltzmann equation (|4.40[) , associated with 
particle production as a consequence of the D-particle foam. This interpretation in terms of particle production is 
phenomenological. As is well known, the concept of particles itself is an ambiguous one in curved space-time. Particle 
production within the framework of kinetic theory can be incorporated as a viscous pressure. It has been shown that 
the inclusion of particle production, in a general relativistic background, leads to a kinetic equation of the form [54] 

p^I, - r" wP v|^ = c[f) + H(x, P ) 

(using standard notation). Here H is the source term representing particle production (or decay) and C collisional 
events. Moreover, in such a framework, it is assumed that C 

• is a local function of / 

• implies number and 4-momentum conservation and 



induces non-negative entropy production and vanishes if and only if / has the form of local equilibrium distri- 
bution 



but otherwise its structure does not need to be specified. The particle number 4-flow, N 1 * is defined as 

N* = j dPp»f(x,p) 

where dP = 5(p v p lJ + m 2 ) — gdp° dp 1 dp 2 dp 3 . H is taken to have the form 

H(x,p)=C(x,p)f°(x,p) 

where f° is the equilibrium distribution; (^(x,p) is parameterized in terms of 4- velocity and momentum as 

£{x,P) = r^r + v(x) 

t{x) 

and two arbitrary functions v{x) and t(x). This form of H is compatible with but more general than the source that 
we derived for D-foam. It can then be shown that [5J| 

N.» = I HdP. 



r = - I HdP 



On writing, 

r = 

n 

(corresponding to our source notation), n being the particle number density, the following relation |54j 

nT = N K + v(x)n 

T 

can be derived. In our case r is large and so can be neglected. Hence a positive v(x) corresponding to particle 
production corresponds to a positive T in (|4.4f I) . It is in this sense that we regard our source T term in (14.41 [) as 
corresponding to particle production. 
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In terms of particle-physics DM models, this implies that in addition to the m 2 x /g x there is another parameter, 
the foam fluctuations cr 2 , which enters the game. In our case, the foam corrections to the effective degrees of freedom, 

— , contribute at most 26% to the final result (14.721) . 

Schematically therefore, we can say that thejaresence of the foam and the requirement of a WIMP-miracle-like 
situation, as implied by the cosmological data 
look like: 



; pre; 
[13], 



fixes combinations of parameters, which in order of magnitude 



9x 



O 



2 <j 
2 9s \ "* A 2 



(4.73) 



for cases where the particle nix of a model, including hidden sector ones, constitutes the dominant DM candidate. 

Whether the D-foam effects on the DM relic abundances are observable depends on the value of the D-particle 
mass scale M s /g s , which in the modern version of string theory is a free parameter. In the next section we proceed 
to discuss briefly the phenomenology of the model, and compare it with some other related interesting cosmological 
models available to date. 



V. COMMENTS ON PHENOMENOLOGY AND COMPARISON WITH OTHER MODELS 



The corrections due to D-foam are in general small, as expected, and they can only be significant for low-string-scale 
models and heavy DM candidates. For traditionally high string scales (M s > O(10 16 GeV)), in order for the D-foam 
effects to be significant one needs superheavy DM, with masses higher than M s /g s . However, the effects of such a 
superheavy DM will be eroded by inflation [28] ; moreover superheavy DM would not be produced significantly during 
a reheating phase of the Universe after its exit from the inflationary period. 

For intermediate string scales (55[, where the quantity M s /g s could be of order fO 11 GeV (which is the order of 
the GZK cutoff of ultra-high energy cosmic rays), there could be significant modifications in the relic abundances 
of superheavy DM particles with masses of this order. Such superheavy DM particles can be produced during 
reheating [56], but in view of our scenario above, their relic abundance will be modified from the Standard Cosmology 
result. The presence of superheavy DM, with increased relic abundances, might provide an explanation for the 
production of at least part of the spectrum of the ultra high energy cosmic rays, with energies of order I0 20 eV. 
Hence, the effects of D-foam on such scenarios are worthy of investigating further, especially in view of the fact that 
the density of D-particles might be significantly higher at earlier eras of the Universe, leading to stronger stochastic 
effects O(Af). 



A. Low (TeV)-string-scale D-particles 

For low string scales, of order a few TeV, the effects of the D-foam on thermal relic densities would be more 
significant. In fact depending on the type of DM considered, the effect could be constrained or falsified already by 
the WMAP data [III, [3], since the induced increase of thermal relic abundance leaves less room for supersymmetry 
in the relevant parameter space. In certain cases, it may exceed the allowed region set by WMAP. Indeed, the DM 
relic abundance determined by the WMAP seven-year-data mean (at 68% C.L.) [14], is: 

(Q^) WMAP,ACDM = 0.1126 ± 0.0036, (5.1) 

thus putting the observational error at the 0(1O -3 ) level. 

In our phenomenological analysis we have found it convenient to parameterize the foam-enhanced abundances (|4.70[) 
as follows: 

O h 2 

"x'^O i,22 

WW) = + (T m * ' 

V X 0/ no source 

„2 3 167T 2 m 2 

a 2 = 1.259 -fe V A 2 , (n x h 2 ) no source = (2.6 x lO" 10 GeV~ 2 ) - 4 x , (5.2) 

where the fudge factor k parameterizes the deviation of the coupling from the weak interaction coupling (/weak that 
characterizes the standard WIMP case [53| and the parameter a 2 encodes the foam fluctuations. 
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The impact of the space-time D-particle foam on the predicted DM relic abundance and the comparison of the latter 
with the latest WMAP measurements are depicted in fig. O We observe that the parameter cr enhances the thermal 
relic density 0, x h 2 , decreasing thus the upper bound of the allowed DM masses compared to Standard Cosmology 
(er = 0). The Vl x ti^ curves are drawn for three values of the factor k, namely k = 0.5, 1, 2, demonstrating the (large) 
dependence on the variation of the coupling constant. 





U.l O 




U.l D 




0.14 


CM 




O 

■C 


0.12 


>< 




a 


0.10 




0.08 




0.06 



_i i i i i i 

k=0.5i 
/ i 

1 7yr WMAP / / 


f « * / II, / 

1 » -* / II* I 
11*1 II* 1 

/; /k=i/ // / 
I*/ ii* i 

/* lit * 

' / / *V / / 


/ i> 
# - * 


/ /// / 7^ 


#**//*/ / 
//*/'// * » 
/// / > * / ' * 

// //// / / -- cr=0.0007 
// /// / s'S —ct-0.0014- 

i MM i M ».* i Mi J> .« 



200 



400 



600 
[GeV] 



800 



1000 



Figure 2: Predicted values of f2 x /i 2 as a function of the DM particle mass m x for various deviations from the weak 
coupling (7woak' k = 0.5 (red), k = 1 (blue) and fc = 2 (green). Different values of foam fluctuations a are drawn: 
cr = (dotted) corresponds to the no-source estimate; a = 0.0007 GeV" 1 (dashed) and a = 0.0014 GeV" 1 (solid) 

describe stochastic D-particle foam models. The yellow horizontal band represents the 7-year WMAP observational 

value (HU) for il x h 2 [H|]. 



In fig. 02 the range for the parameter a has been chosen so that a 2 m x < 0.1, which is a necessary assumption for 
the validity of our perturbative approach on which (|4.70[) and (|5.2p are based. This is evident in fig. [31 where it is 
shown that this condition is met for the whole range of WIMP masses m x . In the same figure, we also present the 
acceptable regions in the (m x ,tr) plane, i.e. those compatible with the 'WMAP corridor" (|5.1|) . for various values of 
the parameter k as defined in (15.21) . It is clear from the plot that there is a shift of the allowed range of m x towards 
lower values and a narrowing of this range with increasing parameter cr. 

For more general models, where the dark matter may even come from hidden sectors and it is unrelated to weak 
interactions, one may use for the DM relic abundance (cf. (|4.73[) and (|5.2[0 : 

2 

n x h 2 oc^f(l + a 2 m 2 x ). (5.3) 
9x 

The effects of foam in such general models are considered in fig. [U where we plot the WMAP-allowed contours [14J of 
ttxh- 2 — 0.1126 in the (mx,gx) plane for various characteristic values of the foam fluctuation parameter cr, including 
the foamless case (a = 0), for comparison. The standard WIMP region, where the coupling g x ~ g' ~ g woa k, is 
also indicated. We observe that, for the allowed values for the foam- fluctuations parameter, < cr < 0.0015, the 
WMAP contour falls within the corridor mx oc g\ , which is a condition sufficient for a model to provide a viable DM 
candidate [53j. 

We note, therefore, that it is a generic feature of the foam-induced Finsler backgrounds to tend to reduce the 
cosmologically allowed range of neutralino (or in general the DM candidate) masses, relative to the foamless cases, for 
fixed low string scales (translated to relatively low — compared to Planck scale — D-particle masses). For low string 
scales (O(10 TeV)), D-foam models may lead to significant modifications of supersymmetry searches at colliders, 
especially in the context of neutralino DM models with the neutralino being higgsino- or wino-like, with masses up 
to 0(1 TeV). 

We also note at this stage that in models where the neutralino DM is higgsino- or wino-like, there may be other 
reasons for an increased relic abundance, for instance slepton co-annihilation [57], and in fact our effects of the foam 
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Figure 3: Regions in the (m x , a) plane allowed by the latest WMAP observations (|5.1j) (lij for various deviations 
from the weak coupling g wca ,k'- k = 0.5 (orange), k — 1 (purple) and k — 2 (green). In the grey area, the foam 
fluctuations are small enough for the approximation (15. 2p to be valid. 




[GeV] 

Figure 4: Contours of Qxh 2 = 0.1126 [1J| in the (mx,gx) plane for various values of foam fluctuations: a = 
(dotted) corresponds to the standard, no-source estimate; a = 0.0007 GeV -1 (dashed) and a = 0.0014 GeV -1 (solid) 
describe stochastic D-particle foam models. The corridor m wea k = {mx / 9x)9' 2 ^ (100 GeV, 1 TeV) (between the 
thin dashed lines) shows the region where mx oc g\. The ellipse marks the WIMP region, where gx ~ g' ~ 9weak- 
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are of comparable strength in some of these cases. Hence, from our point of view, in such models, the addition of a 
low-scale D-foam fluctuating background will lead in general to further constraints on the cosmologically allowed DM 
masses, since the presence of the foam will lead to further enhancement of the DM thermal relic abundances. 



B. Comparison of D-foam with Other Cosmological Models 



Let us now compare the results on the effects of D-foam backgrounds on dark matter abundances obtained in 
the previous sections, with some other, but somewhat related, non-conventional cosmological models. First of all, 
we shall consider the non-equilibrium string cosmologies of [26| . involving time-dependent dilatons 4>(t). The low 
energy effective field theories are (super)gravity theories, with these non-trivial dilaton fields. As discussed in [25j|, 
the presence of the latter, induces source terms in the Boltzmann equation describing the thermal DM relic 

abundances in such Universes, of the form: 

r da (t) = k (5.4) 

where the overdot denotes derivative with respect to the cosmic time in the Einstein frame, i.e. in a frame in which 
the four-dimensional Einstein term in the effective (super)gravity action assumes the canonical form. A self-consistent 
solution of the cosmological equations of these models [26] yields the result < 0, which implies that, in contrast to 
our D-foam effect discussed here, the time-dependent dilaton acts as a sink rather than a source for particles. This was 
important in the foamless cases, since it could lead to a dilution of the DM particle abundance by terms of 0(1/10), 
as compared to the Standard Cosmology models in the absence of dilaton sources, while the baryon density remained 
unscathed. This had profound phenomenological consequences for collider tests of supersvmmetr y I27II . given that 
in that case more room for supersymmetry was left by the WMAP astroparticle constraints [l3], Il4( on DM, and 
thus heavier partner masses are allowed as compared to the Standard Cosmology. Furthermore, other astrophysical 
observations such as type-la supernovae and galaxy ages set tight constraints |58| on the parameters of these models 
through their prediction for the dark energy contribution. 

In the presence of D-foam, on the other hand, cosmologies with time dependent dilatons exhibit competing dilaton 
and foam effects, with the result that, for low string scales, the relevant contributions might cancel each other out. 
For instance, one might obtain a combined time-dependent-dilaton-foam source in the relevant Boltzmann equation 
of the form (cf. (|4.4ip ). which in a schematic way would look like: 

o 2 

r c ombinod ~ + jft \C (Ha\t)mA 2 [9T + 2m]) | , (5.5) 

where A 2 denote dimensionless foam fluctuations. As we have discussed previously, for low string scale the two 
(opposite) sign terms might be of the same order, at least for some range of the DM masses, thereby canceling 
effectively each other out and thus modifying the conclusions of (27| |. Similarly to the dilaton cosmologies, though, 
the foam affects predominantly only the electrically neutral DM particles, and not the charged baryon density, for 
specifically stringy reasons outlined above, associated with either charge conservation (type-IA strings 0, HOy ) or 
compactification details (type-IIB theory |llj). 

Next we come to the models of [30], involving non extensive Tsallis statistics [3l| for particles. As already discussed 
in section ITO1 our D-foam models provide a microscopic origin of effects that mimic in some but not all respects, those 
of a non-extensive statistics a la Tsallis. As in |3(j, the distributions functions that determine the thermodynamics 
and statistical aspects of the Early Universe, differ slightly from the Standard Cosmology ones, by perturbatively 
small terms proportional to the foam fluctuations, assumed small. However, as we have explicitly already mentioned, 
there are major differences from that approach, concerning not only the different functional form of the pertinent 
expressions between the two approaches, but also the fact that the energy appearing in our expressions is not the 
actual energy of the particle but rather its on-shell Minkowski counterpart, as well as, and more importantly, the 
dependence of the corrections on the momentum content of the particle, as a consequence of the Finsler nature of the 
space-time metric deformations, induced by the interaction of the DM particle with the foam defects. 

Another important difference concerns the form of the Boltzmann equation, which in the approach of [3fl | remains 
unmodified. The non-extensive statistics effects can only be seen in the equilibrium expressions of particle number 
densities. Indeed, in terms of the parameter Y, defined above, one has for the Tsallis cosmologies [30] 

dY 

-^L = -(au)s q {Y q 2 -Y q 2 eq ), (5.6) 

where the modified Y q ^ eq can be calculated explicitly using the modified distribution functions. Thus, in the Tsallis 
cosmology cases (30j . the modifications on the DM thermal relic come mainly from the factor |s2 as a result of the 
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g-statistics deformation on the available degrees of freedom, which also affect the freeze-out points of the various 
species, by terms proportional to (q — 1), which in the approach of [30] is assumed perturbatively small. 

In contrast, in our D-foam model, as we have seen above, in addition to a modification in the pertinent particle 
distribution functions, which affect the effective degrees of freedom, the Boltzmann equation is also explicitly affected 
by the D-foam through a time dependent source term T(t) (jl.ip , (14.41(1 . which, because of its positive sign, is expected 
to lead to particle production. In fact, in our case, the foam-induced modifications to the freeze-out point are negligible, 
at least for the weak foam cases we consider here. On the other hand, for WIMP-based models, there is at most a 26% 
correction to the available degrees of freedom from the foam contributions. Thus, the main effects on the DM relic 
population come primarily from the foam-dependent source terms in the Boltzmann equation (|1.1[) . associated with 
particle production and only partly from the corrections to the available degrees of freedom. This is to be contrasted 
to the Tsallis entropy cosmology of [30], where the non-extensive statistics effects are exclusively associated with 
modifications of the effective degrees of freedom of the system. 



C. Other effects of D-foam and complementary constraints 



There are other interesting effects of the D-foam model, which can provide complementary constraints on the 
parameters of the model, by means of entirely different phenomenology. In particular, the non-trivial interactions of 
photons with the D-particle defects in the foam, result in a non-trivial sub-luminal refractive index [lOL |ll|, with the 
anomalous effects being proportional to the first power of the ratio £? bs/-Meff,QG) where -E bs is the (observed) photon 
energy and Megf,QG is the effective quantum gravity scale, given by 

M eff , QG = (5.7) 
n(z) 

with M s the string scale. The function n(z) is the density of D-particle defects in the foam, encountered by the 
photon as it traverses the cosmic distance from emission till observation. This quantity is essentially arbitrary in 
the context of the D-foam models, as it depends on microscopic details of the higher-dimensional bulk physics @, 13 ■ 
Phenomenologically, information on n(z) (mainly constraints on upper bounds) can be provided, for late epochs of 
the Universe, i.e. red-shifts z < O(10), by studies of the arrival times of cosmic photons from distant celestial sources, 
such as Active Galactic Nuclei and/or Gamma Ray Bursts. The idea is simple: there is at present evidence for 
delayed arrival of high e nerg y photons, as compared to their lower-energy counterparts, from a few cosmic sources, 
either AGN [48[ or GRB (49( . The interactions of photons with the D-particles in the foam, cause time delays of the 
more energetic photons by an amount (lol . [llj : 

atob , = r d2 ^2±i), (5. 8) 

Jo M s H(z) 

where H(z) is the Hubble rate, which is a function that depends on the details of the underlying cosmological model. 
By assuming that the foam is the primary source of the delays (which, of course is an assumption, since astrophysical 
effects at the sources could also be responsible for (part of) the delays of the energetic photons), one can obtain 
restrictions on the density of defects and the cosmological parameters of the model, encoded in the functional form of 
H{z). If at late eras of the Universe, z < 10, the Standard Cosmological Cold-Dark Matter model with a cosmological 
constant A > is assumed to be a good approximation of reality, as the current evidence suggests, then one can 
see [59| that fits of all the available data on observed photons delays at present imply constraints on the density of 
D-particle defects n(z), and in fact require that the latter is not constant but drops by almost two orders of magnitude 
as the red-shift climbs from z = 0.03 (MAGIC photon delays 0]), where M s /n(z) = 0(1O 18 ) GeV, to z = 0(1) 
(GRB 090510 observations by Fermi), where M s /n(z = 1) > 1.22M P , with M P = 1.2 x 10 19 GeV the Planck mass. 
In this way the 0(1 min) delays of the TeV photons from the AGN Mk501, observed by MAGIC [H], can be fitted 
with the same D-foam model as the 0(0.5 sec) delays of the 30 GeV photons from GRB 090510 (at z ~ 0.9), observed 
by Fermi. 

This is important information, since it implies that the density of defects in the bulk of the D-foam brane model 
of a a may be non-homogeneous but increases as the time elapses, at least at late eras. Since the density of defects on 
the brane world enters the expressions for the statistical averages over populations of D-particles, which in turn affect 
the induced stochastic fluctuations of the space time, and through them the cosmology of the model [{J [l(|, as we 
have discussed above, one should incorporate such a feature into the analysis. This is left for future work. We stress 
once more, however, that we consider it as an interesting complementary test of the model, providing independent 
constraints on the model parameters, apart from the ones discussed here. 
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Before closing this section we would like to mention another phenomenological consequence of the D-particle foam 
model, namely a qualitative prediction on potential CPT violating modifications [60] of the Einstein-Podolsky-Rosen 
(EPR) correlations of neutral meson pairs produced in meson factories [glj. Specifically, if there is a non-trivial 
density of D-foam defects present during the decay of, say, the $-Meson in Kaon factories, then in general one may 
have a modification of the standard EPR correlators between the long (Kl) and short-lived (Ks) kaon products of 
the decay, which is parametrized by a complex parameter u>, whose amplitude may be estimated to be given by |6l) : 

£ 2 fc 4 
Mj [rag ~ wl) 

where £ 2 is a fudge factor, depending on the density of the D-foam at red-shifts z = (i.e. in the current era) and its 
fluctuations A 2 . However, the fact that neutral Kaons have quark and gluon substructure complicates the situation. 
Indeed, due to their electric charge, the quarks are not supposed to interact dominantly with the foam, only the 
gluons can do this (lOl. [Tl|. This may introduce additional (and significant) suppression factors in £ coming from the 
strong interaction sector of the effective field theory describing the Kaon dynamics. 

The current experimental bounds on uj may be summarized by the results of the KLOE experiment at the Da$NE 
<I>-factory in Frascati (Italy) [62]: 



I cxp 



< 1.0 x 10~ 3 at 95% C.L. (5.10) 



At least one order of magnitude improvement is expected on these bounds with the KLOE-2 experiment at the 
Upgraded DA$NE [63]. Taking into account that \mt — ms\ ~ 3.48 x 10~ 15 GeV, and ignoring possible suppression 
of £ from the strong interactions, we observe from (|5.9I) . that for low string scales M s /g s of order of 10 TeV, for which 
the D-foam contribution to the DM relic abundances (|4.70j) could be significant, the above bound (|5.10|) translates to 
£ < 10~ 12 . For high string scales, where M s /g s is of the order of Planck mass, there are no significant constraints on 
t 

In the context of the microscopic models considered here, one may estimate [64j £ theoretically, for the case where 
there is with certainty one D-particle defect present during the decay of a $ meson, so that the initial entangled 
state of the Kaon products of the decay will be affected by the presence of the foam. In this optimistic scenario, one 
finds |£| 2 ~ (m 2 + m 2 s )/k 2 : where k is the spatial momentum of the Kaon states in the detector. In general, though, 
there are fudge factors in front of this estimate, taking proper account of the probability of a D-particle being present 
during the decay of the initial $ meson, which itself depends on the details of the foam (such as density etc.), as well 
as the non-elementary quark and gluon structure of the Kaon particles, which plays an important suppressing role of 
foam effects, since the foam interactions between the electrically charged quarks and the neutral gluons have different 
strengths, as discussed above [HI EH- Thus, at present, lacking a microscopic non-perturbative detailed estimate 
of these complicated effects, we can only impose bounds on such parameters. Nevertheless, the above diverse tests 
indicate how in general space-time foam models of the type considered here may be constrained in the near future. 



D. On the production of TeV- mass-scale D-particles at Colliders 



We would like now to comment briefly on the production of D-particle defects in particle collisions at high energy 
colliders, such as the LHC. In general, the production mechanisms will depend on the type of D-particles (and on the 
associated microscopic string theory model) considered. In this section we would like to make the discussion more 
general. 

We will first distinguish our D-particle cosmological model from the generic D-matter analysed in [2(j. In the 
latter work, D-particles are viewed as ordinary effectively point-like excitations of the vacuum. They arise from the 
appropriate compactification of Dp' branes around p' cycles, with radii of the order of the string length £ s = 1/M S . 
Since they have non-trivial couplings with Standard Model particles, they behave like standard dark matter candidates; 
hence their production mechanism at colliders will not differ from those of conventional dark matter. In fact in [20], 
such D-particles have been treated as providing a leading candidate for dark matter in our Universe. As in our case 
above, the mass of these D-particles is of order M s /g s , i.e. TeV scale defects can be present in low scale string theory 
models. 

In general brane defects have spin structures; hence stable D-matter states can be either bosonic or fermionic, 
corresponding to the bosonic or fermionic zero modes of stable D-branes respectively. In our cosmological model 
described in previous sections, we treated the D-particles as background defects, where their spin was not relevant 
for our induced effects; on the other hand, when we consider issues like the production of localised branes, then their 
spin is relevant, since the corresponding interaction cross sections, for instance with Standard Model particles, such 
as nucleons, do depend on spin. In this respect, the production of localised D-branes by the high energy collision of 
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Standard Model particles at colliders can incorporate even more complicated charged defects, consisting of various 
configurations of compactified Dp'-branes around p'-dimensional manifolds. 

In a low-energy, string-inspired, effective field-theory action, the leading interactions of the D-particles with Standard 
Model matter are provided by terms with the generic structure (omitting Lorentz derivative or Dirac-matrix structures 
for brevity) [H 



oc gu D D (gauge bosons) 



(5-11) 



The symbol oc in front of each type of interaction is included to denote form factors that arise from tree-level string 
amplitude calculations with appropriate Dirichlet and Neumann boundary conditions for the open strings [20]. The 
relevant string amplitudes have three- vertex insertions on world-sheet discs, with the appropriate boundary conditions: 
Dirichlet for the vertex operators describing an excitation of D-particles, and Neumann for ordinary open strings 
describing (ordinary) gauge bosons. The resulting effective three-point amplitude is expanded in powers of the ^fg~ a 
of the open string coupling g s < 1 , assumed perturbative, 

g 1 J 2 F{s,t,a , ) + g s g(s,t,a , ) + ... 

1 /2 

where s,t are Mandelstam variables. In the context of string/brane theory, the Yang-Mills coupling gym ~ 9s ■ 
Thus, the effective coupling gu is proportional to gym, but renormalised by appropriate kinematic factors and higher- 
order corrections in g s . For momenta of the same order as the low-string scale the corrections may be important 
at high energy scattering. Unfortunately string corrections of D-branes are not understood very well at present to 
be able to give any quantitative prediction of the effective coupling go, especially for low scale D-particles; for our 
purposes we can only use generic effective field theory arguments in order to study the phenomenological features of 
the D-particles. These are in general dependent on the microscopic-model. The terminology "gauge bosons" in (|5.1 1[) 
will represent not only gauge bosons but also other Standard Model particles. As a result of (|5.11[) . for instance, 
one may have the graphs of fig. I5al arising from quark/antiquark scattering. The D-matter/antimatter pairs can be 
produced by the decay of intermediate off-shell Z-bosons of the Standard Model, which is in agreement with f)5 . . 
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Figure 5: (a) Feynman diagrams at the parton level for the production of D-particles by, say, qq collisions in a 
generic D- matter low energy model of [2(j, based on (I5.11[) . This is only an example. There are many other processes 
for the production of D-matter from standard boson decays or gauge boson fusion, which we do not consider in our 
qualitative discussion here, (b) Production of conventional dark-matter particle-antiparticle (xx) m effective field 
theories, assuming that dark matter, which may co-exist with D-matter, couples to quarks via higher-dimensional 
contact interactions [65, 66]. (c), (d) The dominant background processes within the Standard Model framework. 



The D-matter pairs produced in a hadron collider will traverse the detector and exit undetected as they are weakly 
interacting only; however there will be large transverse missing energy (E^ lss ). Hence analyses requiring high Elf ms 
and one energetic jet would be of high relevance. The dominant background from Standard Model processes involves 
the decay of Z to a neutrino pair and of W + to a lost ("invisible" ^j nv ) lepton and a neutrino (cf. figs. l5cl and [5dl 
respectively) . Such searches have been performed by CDF at the Tevatron |67l | and by the ATLAS [68] and CMS |<39j ] 
experiments of the LHC at centre-of-mass energies of 7 TeV giving null results. The interpretation of the such analyses 
may be deployed to constrain D-matter parameters, as will be discussed below in the context of our specific D-foam 
model. 
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In the approach of [20|, the D-matter behaves as the dominant species of dark matter in the Universe; they do 
not consider (as we do here) the effect of D-matter in the bulk region of the brane world. The cross section for the 
above production of DD-pairs, a(qq — > DD + ...), is related to the total Dark-Matter annihilation cross section 
ftotai annih into Standard Model particles via the principle of detailed balance (which one may assume to characterise 
string theory effective action models), DD -> //, where / denotes Standard Model particles (including quarks of 
course). The latter cross section can be constrained by cosmological data, if the D-matter is assumed to play the role 
of the dominant species in the Universe, since the total annihilation cross section is inversely proportional to the relic 
abundance of D-matter, floh 2 oc j— - j 1 — p^-, in a standard notation. It is for this reason that the mass parameters 
and couplings of the D-matter have to be constrained in such a way that the region of the cross section is near the 
'WIMP-miracle' region. In this approach one also considers the scattering of the D-matter against nucleons, which 
provides, as usual, complementary constraints on the dark matter parameters (especially its mass and couplings) via 
direct searches. Needless to say that the relevant cross sections depend crucially on whether the D-matter is bosonic 
or fermionic (D-particles can have integer or half-integer spins, as a result of the corresponding bosonic or fermionic 
zero modes they characterise D-matter). 

However, it is possible that, the D-matter particles being relatively heavy and thus slow moving, deposit all their 
energy inside the detector. Provided the deposited energy from the D-D-particle pair is large enough, this may lead , 
for instance, to other observable effects in the MoEDAL detector of the LHC collider [7(j, which is a plastic nuclear 
tracks detector designed to search for magnetic monopoles and other highly ionising particles that may characterise 
several models of new physics, such as supersymmetry. Despite being electrically neutral, D-particles can operate 
in a similar way as a monopole, by destroying the chemical bonds in the plastic sheets that surround the MoEDAL 
detector, thereby leading literally to the formation of holes. In fact D-particles are analogous to 't Hooft-Polyakov 
monopoles (solitons) in string theory. One difference is that the couplings of D-particles may be perturbative, in 
contrast to the monopole case, which has a coupling oc 1/gyM- From the shape of the various holes along the path 
inside the plastic sheets one can in principle determine the mass of the heavy particle and its type. 

Finally a third way of producing TeV-mass D-particles at particle colliders would be through the formation of TeV- 
scale black holes, which have short life times and then decay, a la Hawking, giving rise to pairs of particle-antiparticles 
that describe localised excitations of the low-energy effective field theory along with stable D-particle-D-antiparticle 
pairs. 

As we have indicated, our model for D-foam has some essential differences from the generic D-matter models of [20]. 
First, in our type IIA models of string foam |8l-ll0| the D-particles are truly point-like and the D-foam is considered as 
a stochastically fluctuating background, with no (or suppressed) couplings between the D-particle defect and Standard 
Model fields. For us, the dominant interactions of D-foam with ordinary matter are gravitational in nature, provided 
by the distortion of space-time due to the recoil of the defect during the capture/string-splitting processes. In this 
sense, D-particle in our models cannot play the role of dark matter in the Universe. Their presence and gravitational 
interactions (through recoil) with "conventional" dark matter candidates, affect the relic abundance of the latter. 

In these type of models, D-particle defects can only be produced at colliders by first forming TeV-size black holes in 
low scale string models, which then undergo Hawking radiation, in which pairs of TeV-mass D-particle/anti-particle 
are produced along with all other Standard Model or other exotic particle pairs, depending on the model. In this 
sense, for this type of string IIA D-foam [H, [l(|, the only collider signatures are the indirect ones via the increased 
relic abundances of dark matter, as discussed in the previous section. 

In general, however, the D-particles (arising from compactification) in other types of D-foam models, such as 
type IIB string theory models of [ll|, can couple to the Standard Model fields, but with couplings suppressed in 
comparison with the standard model ones. In such models one may consider the formation of D-particles in some 
high energy collisions and in the early universe, along the lines of D-matter discussed previously [20]. Nevertheless 
in our D-foam cosmology, as we have emphasized, in contrast to the models of [20], the D-particle contributions to 
the Universe vacuum energy density can be compensated by interactions of D-particles in the bulk regions of our 
configurations (cf. fig. [IJ. In section PV A[ it was noted that there are negative energy contributions from bulk D- 
particle defects that cancel out (on average) the positive energy contributions due to the D-particle masses and relic 
densities (cf. I|4.7p ). Hence, in our models of D-foam, in contrast to the generic D-matter considerations of ref. (20| . 
there is no restriction in the relevant annihilation cross sections into Standard Model particle/antiparticle pairs, //, 
<7 a (DD — > ff). Such restrictions apply only to the conventional dark matter annihilation cross sections that do exist 
in our D-foam cosmology. 

Let us give a concrete example of D-particle production at colliders in our type IIB string models of D-foam, 
considered in ref. (TT| . Here D-particles are essentially D3-branes wrapped up around appropriately compactified 
three cycles (e.g. tori) [ll|. Such defects will look "effectively" point-like from a four-dimensional space time point of 
view, but the difference from the truly point-like defects of the models of pj-Ho| lies on the fact that in this case there 
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are non-trivial (dimensionless) couplings <?| 7 between the D-"particles" and the Standard Model ones [Tl| 

a (1-55) 4 , 2 (1.554) 4 2 _ 2 ,_ 19 . 

where V43 is the average three-spatial-dimensional volume which contains one D-particle (assuming uniform distribu- 
tion of defects in the foam) , R' is the compactification radius of the dimension transverse to the wrapped- up-D3-brane 
(D-"particle") and gj denotes the (dimensionful) gauge couplings of open strings leaving in the D7-brane in the model 
of which after appropriate compactification describes our observable world (the corresponding dimensionless 
couplings are defined as 177 = g?/£*). Hence, gr is proportional to the Yang-Mills couplings (gym oc i/g~7, where g s 
is the (dimensionless) string coupling) of the various Standard Model gauge bosons, say. It is important to notice 
that the form (|5 . 1 2|) characterises exclusively the D-foam model of [H| . and should not be confused with the generic 
models of D- matter discussed in [2(| . As discussed in [ll|, in order to avoid observable Lorentz-violating effects in 
the dispersion relation of electrons that would affect the synchrotron radiation spectrum from distant nebulae, such 
as Crab, by effects that are not observed experimentally, one must have 

Such a bound, for instance, is easily satisfied in our D-foam string models with Va3 ~ (10^ s ) 3 , R' ~ 3.4 x 10 2 £ s . 

In view of ()5. 12[) . and assuming, for definiteness that D-particles are fermionic, one has trilinear interactions between 
the type IIB-D-"particles" with ordinary Standard Model particles, such as gauge bosons, of the form (|5.11|) (omitting 
for brevity Lorentz derivative or Dirac-matrix structures for brevity) Ho| 

oc 537 D D (gauge bosons) , (5-14) 

where, as in (15 . lip , the symbol cx in front of each type of interactions denote kimematical form factors that arise from 
tree-level string amplitude calculations. In view of (|5 . 13|) , in this model, the trilinear coupling of D-particle/ antiparticle 
pairs with the Standard Model gauge bosons would be much weaker than any typical Standard Model coupling. Indeed, 
the mass scale of the effective four-fermion (q qD D) interaction, representing the mediating off-shell Z-boson entering 
the D-particle production process in fig. [Sal nas the form 

M«=^=~M*TrV* (5.15) 

V53757 97 

where Mz ~ 90 GeV is the mass of the Z-boson and gj represents the Standard Model weak interaction coupling on 
the brane world. Thus, on account of (|5.13[) . M e g can be of order of a few TeV, making the D-production graph of 
fig- EI] °f similar order to the conventional dark-matter production graphs in the model (c/. fig. I5b|) . However, in our 
case, unlike standard scenarios involving only a single type of dark matter dominant in a given channel [65], |66j, we 
have a mixture of dark- matter contributions, coming from D-matter and conventional (WIMP- type) dark matter (%) . 
As a result, the detailed phenomenology is more complicated, and depends on the relative masses of stable D-matter 
and conventional dark matter. In the case of low-string scale (TeV) D- and dark- matter, one is in the border-line 
region for defining the low-energy effective low-energy string action, and stringy corrections are important in our 
model. 

The null searches of, say, ATLAS detector in finding missing energy + jets, will thus impose, in the context of our 
model, a bound on the dark matter mass m x and D-particle couplings (737 (and thus 77, and therefore the D-foam 
density). With the available LHC data at y/s = 7 TeV and 1 fb _1 integrated luminosity, the current lower bounds 
on the mass scale below which the effective field theory of dark matter is valid, is 0(200 GeV) [66]. Therefore we 
need the full LHC potential (1 fb _1 at ^/s = 14 TeV) to produce, and thus constrain the parameters of, our TeV-scale 
D-particles (c.f. ([535]) ). 

It should be remarked at this stage that, since r\ depends on the density of foam, the larger the elementary ("unit") 3- 
volume containing a D-particle or the radius i?', i.e. the smaller the D-foam density, the smaller the effective coupling, 
and hence the weaker the relevant interactions. In this case, the D-particle interactions with the Standard Model 
excitations can be purely perturbative. Unlike monopoles, which are non-perturbative for weakly coupled strings, 
the D-particles can be treated perturbatively. In fact, this similarity of the D-particles with the 't Hooft-Polyakov 
monopole, makes the MoEDAL detector of the LHC [70] a natural place to look for their potential signatures through 
post-production distinctive signatures ("etching") on the plastic sheets of the detector. In this way, the non observation 
of such D-particle pairs (with masses in the TeV range) at high energy collisions can lead to bounds on their couplings 
(I5.12p . and on the density of D-particles. 
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It must be stressed that in our cosmological models of D-foam, the D-particle defects do not constitute the bulk 
of the dark matter. Thus, their annihilation cross sections, which are proportional to (<?37<77)/M|, = (.937(7733 )/M^ : 
do not have to be in the 'WIMP- miracle' regions, and can be much smaller for TeV-mass D-particles, Mo ~ M s /g s 
(corresponding to compactification radii of order of the string length). Indeed, as we discussed in section ["IV Al 
as a result of the existence of D-particles in the bulk regions of our configurations (c/. fig. [T|), there are negative 
contributions to the energy density of our brane world, induced by such bulk defects, which cancel out any positive 
energy contributions to the Universe's energy density coming from frozen D-particle relics (c/. (|4.7p V Hence there is 
no restriction in the relevant annihilation cross sections <r a (DD — ¥ //), where / are Standard Model particles, say. 
Such restrictions apply only to the conventional dark matter annihilation. 



VI. CONCLUSIONS AND OUTLOOK 

In the present article we have computed Dark Matter thermal relics (via the appropriately modified Boltzmann 
equation) in the presence of a particular type of stringy space-time foam, including brany D-particle defects as the 
foam structures. We have seen that, as a result of the interactions of the DM particle(s) with the defects in the foam, 
there are induced Finsler-type metric distortions, depending — in addition to the space-time coordinates — also on the 
pertinent momentum transfer. These momentum-dependent metrics have profound effects on the geodesic equation 
of the DM particles, leading to modifications in the pertinent Boltzmann equation determining the (thermal) relic 
abundances. 

The foam acts as a source of particle production [7|, but also affects particle statistics in a way reminiscent of 
non-extensive statistics of Tsallis type. In this latter respect, our D-particle foam models may be viewed as providers 
of microscopic models of non-extensive statistics in Cosmology. The net result, is that the DM relic abundance in 
the presence of the foam is larger than the corresponding one in Standard foam- less Cosmology (ACDM model). 
The effects of the foam in such DM relic populations depend on the value of the string scale and coupling. For 
conventionally high string scales, close to 10 18 GeV, the effects are negligible. For low string scales, on the other hand, 
the foam effects may be constrained by WMAP data on colliders, as they imply considerably less available room in 
the parameter space for supersymmetric particle physics models than in the ACDM scenarios. 

In this latter respect, the D-foam acts as a competing effect to those of time dependent dilatons, characterizing 
certain non-equilibrium stringy cosmologies. Indeed, in such models, the time dependent dilaton fields result in 
negative sources (sinks) in the Boltzmann equation determining the DM relic abundances. In the absence of foam 
effects, this may lead to a dilution of DM abundances by as much as 1/10 of those in the Standard Cosmological Model. 
This would have profound effects on collider, and in particular LHC, phenomenology. The presence of low-string scale 
D-foam, on the other hand, acts in the opposite sense, since it tends to increase the relic abundance. Thus, in the 
presence of D-foam, the time dependent dilaton effects may be neutralized in some models. 

Complementary constraints on the D-particle foam may be provided by high energy gamma-ray astronomy, via the 
induced vacuum refractive index that such models entail. Moreover, interesting tests can be made at particle interfer- 
ometer experiments, such as neutral Kaon or other neutral meson factories, as a result of the modified entanglement 
properties of neutral mesons, due to their non-trivial interaction with the D-foam. 

Production of TeV-mass D-particles at high-energy colliders, such as the LHC, has also been briefly discussed with 
the conclusion that such a production could be possible and lead to observable signatures characterised by missing 
transverse energy only when the collider operates at its full energy. On the other hand, such relatively light defects 
when produced inside the MoEDAL detector of the LHC may deposit sizable amounts of energy via collisions with 
the molecules of the nuclear track plastic arrays, thus being detectable through the induced tracks. 

Needless to say there is still a lot to be done before definite conclusions are reached on the falsification of models 
dealing with the quantum structure of space-time. Nevertheless, certain models, such as the D-foam examined here, 
may be falsified or constrained significantly by a plethora of diverse cosmological or astroparticle physics tests in the 
foreseeable future. 
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APPENDIX A: Distribution Functions in non-Expanding Foam Universes 

The number density n and energy density p of species in our stochastic space-time foam, when one ignores the 
Universe's expansion, are given by the basic formulae: 
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((n))d 3 k (6.1) 



and 



P = j£f j ((™r))d 3 K (6.2) 

with ((n)) given in (|3.19p in the text. 

The evaluation of the momentum integrals involved is facilitated by going to spherical coordinates, i.e. setting 
ki = fcsin#cos0, k 2 = fcsin#sin</> and k% = kcos9, where k = \k\ = y/k 2 + fcf + fcf and 9 and <f> are the ordinary 
polar and azimuth angles running between [0,tt] and [0, 27r] respectively. Then, we obtain d 3 k = dk k 2 sin 9 d9 d(f> and 
therefore the angular part of the integration of the different powers of kj appearing in (|3.16p yields: 

d6sm6 = 4ir, (6.3) 
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We commence our analysis with the number density calculation. The integral in (16. ip can be written in the following 
compact form: 

((n))d 3 k = (4tt)J(0, l,m,g) 

y /? 3 (P(l, 1, m, - £P(1, 2, m, 0) + y /? 4 (.7(2, 1, m, - 3£J(2, 2, m, + 2£ 2 J(2, 3, m, £)) } , 



where we have set for brevity: a 2 = (cr^ + a 2 2 + ct 2 ^) with a^i — ^ (we note that the final results will be expressed 



in terms of the parameters aoi) ail d defined: 
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f°° E 2 (E 2 - m 2 ) l+1/2 , 
P(l,n,m,0= dE \ L— . 6.7 

The above formulae (|6.6l6.7p are used to calculate the number densities of the various particle species, either 
relativistic or not, at different eras of the evolution of the Universe. 

In the relativistic limit (fc^T 3> m), we can consider m — > 0, and therefore the integrals (|6.7p can be rewritten as: 
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Hence: 

((n))d 3 k = (4^)/3- 3 / 2 (l;M/3,C) + 



\ y h/5(l; MA + e/s(2; 0] + y [/e(l; - 3£/ 6 (2; ^8, £) + 2£ 2 / 6 (3; ^, £)] | , 



where the functions f s (1; £), f s (2; /J./3, £) and / s (3; /z/3, £) have the explicit form 



we find 



and 



(6.9) 



/ 5 (l;x,l) = -r(s + l)Lii +5 (-e x ), 

/ s (l; a:, -1) = r(a + l)Li 1+s (e x ) , (6.10) 
with Li n (z) the polylogarithm function: Li n (z) = X^fcLi defined for |z| < 1. From the recursion formula 

f a (n+l;x,£) = \f s {n;x,£) - -7-7- is (n;x,£) , (6.11) 



/ fl (2; = (a + 1) (Lt 1+S - Li s (-&*)) (6.12) 



/ s (3; x, = (a + 1) (^Li s - L* 1+s - X -Li s ^ . (6.13) 

In the ultra-relativistic limit /i <C fee? 1 , and in leading order we can set ///? = in (|6.9p . yielding finally the 
relativistic number density for bosons and fermions: 



n b = |^r(3)C(3) + ^ T 5 a 2 C(5) (-y r(6) + y r(7 . ) ((,.14) 

and 

»/ = ||^ r(3)c(3) + 8^ t5 ^ {y r(6) [_27?(6) +7?(5)] + f r(7) [6??i7) (),/,<,i + ^ • 

where £ and are the Riemann zeta and Dirichlet eta functions respectively and we have set: (Tq = cr^. 

The first term in both (|6.14[) and (|6.15l) represents the standard result whereas the second term gives the D- 
particles foam correction. This correction, does not have the usual scaling with temperature T 3 but scales as T 5 
instead. However this is not really in contradiction with Standard Cosmology since the correction, weighted by the 
small parameters croi, can still be far below current observation limits. 

In the non-relativistic, large-mass limit, the quantities £ can be neglected in the denominators of the integrands in 
(|6.7p . that now reduce to the integrals (upon integration over x = (3E): 

g(l,n,m,0) = J dxx (x 2 - (m/3) 2 ) + (exp {x - m/3))~™ 

1 • >)„;2. l (^ 



021+3 
1 



t^-m 2+l ( -) K 2+l (mnf3)(3 2+l e^ n , 
ir \n J 



h(l,n,m,p) = ^+4 J dxx 2 (x 2 ~(mf3) 2 y^ /z (exp (x -up)) 

i+i 



_L_(^y + ' Jri+ ,(mnfflI^li4(l + ^.*» 



where K v is the modified Bessel function of the second kind of order v. 



(6.16) 
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In terms of these functions we obtain for the heavy-species number density p. 211) : 

= gfg r 3 {(47r) 3 (0, 1, m, /3)e^ + ^P [g(2, 1, m, /3)e^ - 2fr(2, 2, m, /?)e 2 ^ + £ 2 5 (2, 3, m, /3)e 3 ^] 



y ~al [h(l, 1, m, /3)e^ - 2, m, /3)e 2 ^] - ^-P [g(2, 1, m, fle"" - &(2, 2, m, /3)e 2 ^] 



where, in the subleading contributions, £ = +1 applies to fermions and £ = — 1 to bosons. In deriving (|6 . 1 T[) we have 
used the asymptotic limit of the Bessel function for large x values, K v (x) ~ \/5 e ~ X i a consideration valid in our 
non-relativistic approach (x = m/3 3> 1). 

In order to obtain the energy densities (|3.22[) in this model, one needs to calculate the average over the statistics of 
the product of the distribution function n r as given in (|3~T6ll and the energy lu t p. lip . By keeping again powers up 
to r 2 and dropping riTj terms for i ^ j, we find: 



3 



exp(2^Cg - m)) 
(expC8(£ - m)) + : 



AfE^tfr*. (6.18) 



Averaging over rj and integrating over the angular part (in spherical coordinates) of d 3 k , yields: 
((n r uj r ))d 3 k = / 87rdfcfc 2 ^-- 



o t2exp( i 9(S-M))+^ 

x 2q2/ 1 k 2 E (3 expQ3(£ - M ))fc 4 exp(/3(£-^)) 2 2 
_cr oP To TTTTTJTz? TvTTT ¥ 7IZZ7~E7~B TvTTTTa To" 71 7^7~S rTTTTTa ^ 



, 12 exp(/3(i? - M )) + £ 5 (exp(/3(£ - /x)) + £) 2 12 (exp(/3(£ - /x)) + £) 2 
/3 2 exp(/3(E - m)) ■ 2 exp(2/3(g - A 1 

"5 (exp(/3(i?- M ))+0 2 5^ (ex P (/3(£- M ))+0 3 7J ' 1 j 



The reader should recall that a 2 = (3 2 ct^ 



In terms of the parameter E 2 = k 2 + m 2 , (|6.19p takes the form: 



((n r oj r ))d 3 k = 8tt(3- 4 \ -P(0,l,m,0 



-ip(2, 1, to, + i^P(2, 2, m, + |P(2, 1, to, - g£P(2, 2, m, + j^ 2 P(2, 3, to, 



(6.20) 



where we have set: 



x (x 2 — m 2 r2 5 ' ' 



j(3 2 y 

J(l,n,m,£)= I dx-^— n , (6.21) 

\p {exp{x - fi/3)) + £) 
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r°° x 2 (x 2 -m 2 B 2 Y +1/2 
P(l,n,m,0= I dx-±- j:/ -» (6.22) 



and 

roo x 3 ( x 2 _ m 2 ( g2\ i + 1 /2 

Q(l,n,m,0= dx- n = J (I + 1, n, m, + m 2 J(l, n, m, £_). (6.23) 

J mf3 (exp(x - nP) + £) 

The relations between J(l,n,m,£), P(l,n,m,(,) and the corresponding quantities J(l,n,m,£) and P(£, n, m, £) 
defined in (|6.7[) are: 

J(l,n,m,£) = J(i,n,m,£), 

P(l, n, m, = n, m, 0- (6.24) 

Now from (|6.20p one can calculate the pressure and energy density for any kind of species at any era of the evolution 
of the universe by calculating appropriately the integrals (|6.2ip , (|6.22p and (|6.23p . 

The non-relativistic pressure density p, for a heavy dark matter species of mass m 3> \k\ = k, is obtained from the 
general expression for the pressure 

^(2^/0^ (6 ' 25) 
upon taking the limit m 3> k, i.e. uj r ~ m to leading order, and using as definition of temperature 

( .,_ .,_ j n k 2 d 3 k nmT. (6.26) 

(the reader is invited to compare this with the definition in eq. (|4. 38[) in the text, for the case of an expanding 
Universe). This implies for the non-relativistic pressure p = + . . . , where the dots denote foam corrections of 
order a 2 , suppressed also by kj/m <C 1 terms. 

Similarly, the non-relativistic energy density is found to be: p — mn + 0(a 2 ). For isotropic foam situations, the 
(small) foam corrections can be easily seen to be cast in the form p = mn(l + Tma^), where the above-mentioned 
temperature definition (|6.26p has been used. 

From these considerations, for cold dark matter species T <C rn, the pertinent equation of (approximate) dust 
emerges: w — p/ p — (l + O(crJ)) <C 1 to leading order in the small foam corrections. 

In the case of relativistic matter, the energy density is found to be: 

P = ^ 1^-^/3(1,^,0 + ^" 6 (-^/ 5 (1,/3M) - \h{^M + |e/a(2,j9A*,0 - -^/e(l,M 

+l£/ 6 (2,/3 M ,0 - J/r(l,/3/i,0 + \Ur%fat) + ^/r(l,M - \Ur%fat) + ^ 2 /t(3, fa 

(6.27) 

Calculating the pressure of a relativistic gas, using the general formula (|6.25p . and comparing it with the expressions 
for the energy density (|6.27p above, it is straightforward to observe that the presence of the foam leads to deviations 
from the standard equation of state for radiation u> ra d = ^, by terms proportional to the foam fluctuations 0(<j 2 ). 
This is an important difference from the case of the non-extensive (Tsallis) statistics cosmology of [30j, where the 
relativistic equation of state retains its conventional form w^ 11 ' 8 = 1/3. 

Upon using relations (|6.10l) . (|6.12p and (|6.13p and setting: p = 0, we find from (|6.27l) : 

Pb = |^r(4)C(4)T b 4 + 9 ±o 2 Tt (^r(6)C(6) - ^r(7)C(6) + ir(8)C(6)) = ^g,n + Hb^dffi, (6.28) 
for bosons and 

Pf = |^r(4)ry(4)T / 4 + 2L<%Tf ^T(6) V (6) - ^(7)^(7) + ^r(7)r?(6) + ^(8)^(8) - jjr(%(7) + ^(8)^(6)) 
f 7 \ n2 ^4 , 91-2^6 ^ 18941tt 6 50841 127i 8 \ /7\ tt 2 4 793.32 _ 2m6 
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for fermions. These equilibrium distributions will be important in estimating the correction to the effective number 
of degrees of freedom. 

Changes to distribution functions can have cosmological implications in principle and we considered them for the 
freeze out of dark matter candidates. 

APPENDIX B: Geodesies Equations in Stochastic Finsler Space-Times 

The geodesies are obtained by using the standard technique of the Euler-Lagrange equations for the action S given 

by 

f T2 pr . , , f dx a \ dx^ dx V 

b = I \f Ldr with L 

leading to: 



' dr ) dr dr 



d dL dL 

= 0, (6.30) 



since r is an affinc parameter. This yields: 



d 2 x» dx^dx" 1 dx^dx^dx? 

N ^ = - T ^^~ 2 9 9,iU:a ' p ~dV~dT~dV (6 ' 31) 



with the convention g^ v . a = -g^g^v and g^. a = gfs-flV- Here, 

.dx v , 1 Ka dx^ dx v 

Equivalently. equation (|6.31j) is written as: 



N£ = s; + g Ka — (g^+g^-.a) + ^g Ka g^. a ^— — . (6.32) 



d2 x X (N -ix\^ n dx» dx v 1 „„ dx^ dx v dx 13 

H^ 2 



We should stress that this is a general modification of the geodesic equation for an arbitrary dependence of the metric 
on v. In our case however, things simplify a lot since g^-.p is zero in all cases except for terms g^i-i = <7j0:i = mri<i 4 (t) 
and also gp V :a\n = 0. The inverse metric g^ v is readily calculated and is given by 

oo a 2 (t) Ql k iri a 2 (t) u l + a 2 (t) k 2 r 2 k^r^it) 
9 = , 9 = , 9 = , 9 = , 

q q q q 

with q = a 2 (t) [a 2 (t)£i*frf + 

We need (TV -1 )^. For a general metric, this can be written as: 

{N-% = 6Z-UZ+(U 3 )l, (6.34) 

where 

dr" 1 drfi dr" 

jjk = g ^ax_ {g ^ a + ^ + ^ gKag0ua ^ax_^_ (6 35) 

This expression may seem complicated but we remember that in the end we always work up to (and include) terms 
0(r 2 ). The other ingredient in the calculation is the non-zero Christoffel symbols which are found to be: 

r o _ 2a(t)a(t) [E t r?k 2 ] _ _ r t fc t a(f)q(t) _ a(t)a(t) 



1 + « 2 W E, r^l} ' ° l l ° 1 + « 2 W Ei r?*?] ' " 1 + « 2 (*) [Ei r 2 fc 2 ] ' 
. _ 2Hr i k l , t _ a{t)a{t)r 2 k 2 t r^rj fcj-a(t)d(t) 
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• rikia{t)a{t) t rik t a{t)ii{t) , rjkja(t)a(t) 



where i? = is the Hubble parameter. We also note that in what follows, the use of summation convention is 
generally eschewed and on the occasions that it is used, it will be clearly indicated. 

After detailed calculations, where only terms up to order r 2 have been kept and cross terms of the form r^j for 
i =/= j have been dropped, equation (|6.33[) . for fj, = i = 1, 2, 3, simplifies to: 

2 .... „ 2 



-Hp'p - —a{t)a{t)r lP * (p°Y + —a\t)a{t)rl (f)p° 

d,T z m z m z v ' ni 



+ ^a 3 (t)a(t)np* J2 {p i f + -^a 3 (t)a(t)r^ (p ) 3 - ± a 5 (t)a(t)r^p° ]T {p>f . (6.36) 

TYlj f~f~L 1~t~L 

3 

These equations are used in the modified Liouville equation that determines species abundances in the stochastic 
D-foam space-time. 
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